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ABSTRACT 


In  recent  years  the  public  has  becono  acutely  aware  of 
the  energy  crisis  that  has  been  inevitable  since  man  first 
became  dependent  on  finite  fossil  fuels.  Today  there  is  emphasis 
on  developing  "alternate"  energy  sources.  The  use  of  v/ind  pov/er 
, for  irrigation  of  farm  fields  is  investigated  in  this  paper. 

Many  options  have  been  proposed  for  harnessing  wind  power 
for  irrigation,  /lore,  two  options  are  investigated  using  a 
system  dynamics  approach.  The  first  option  uses  mechanical 
enei'gy  extracted  from  the  wind  as  input  to  a mechanical  pump. 

A flyv/fioel  servos  ns  a buffer  and  as  storag'o  for  low  v/ind 
periods.  The  second  option  user  the  mechanical  output  from  the 
v/i rich,  ill  to  generate  electricity.  This  electricity  is  then  fed 
into  an  electrolysis  cell  to  yiroduce  hydrogen  gas.  This  gas  is 
then  burned  in  an  internal  combustion  engine  which  drives  the- 
pun-p,  Surplus  hydrogen  is  stored  in  depicted  natural  gas  wells 
for  ure  during  low  wind  periods,  A back  up  system  is  provided 
in  both  options  for  when  the  wind  is  not  blowing  and  the  stored 
encr-gy  has  been  exhausted, 

1.965  v/ind  data  fi‘oi:i  tlie  U.Sf  V.'eather  Bureau  Office  in 
Dodge  City,  Kansas,  is  used  i.n  tiio  simulation  of  the  two  proposed 
systems,  1965  wind  data  is  used  because  it  tends  tc  be  a "typical" 
wind  year.  The  summer  montlis,  duri^rig  which  the  irrigation  will  be 
taking  place,  have  average  wind  velocities  within  one  half  a 
standard  cinviation  of  a 25  year  average  for  the  area.  The  other 
nine  months  each  have  average  witid  velocitior  v/ithin  tv/o  standard 
devir lions  of  a 25  year  average. 

First,  the  general  concept::  of  system  dy;i.'mics  arc  explained. 
Th(n,  models  of  tlic  two  options  are  formal at('d  and  tested  by 


simulating  their  performance  over  a "typical"  year  using  the 
DYNAMO  compiler.  The  results  of  these  simulations  are  then 


compared  with  each  other  and  with  the  current  methods  of  farm 
irrigation. 

Finally,  a glance  to  the  future  is  proj.'osed  which  reveals 
that,  though  the  systems  proposed  are  presently  more  expensive 
than  current  forms  of  energy,  tliis  very  well  might  not  be  the 
case  in  10  to  15  years. 
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CllAWER  I - INTROIHICTION 

The  amount  of  land  ai-ea  bcitvi  firmed  per  person  in  the 
United  States  has  been  on  tlic  decline  for  the  pant  several  years. 
Yet  today  the  people  of  America  are  eating  better  than  ever 
before.  This  is  duo  in  part  to  the  development  of  now  strains 
of  crops  that  reach  maturity  factor  and  £;ive  a higher  yield  than 
the  old  strains.  But  at  the  came  time  the  amount  of  enor»;y 
required  per  p:iven  area  of  farm  land  has  inci'oascd.  Lar(;er,  more 
powerful  pieces  of  farm  equipment  are  beiivi  used  today  than  over 
before,  I t is  customarily  the  case  that  these  machines  are 
fueled  by  an  oil  product.  It  is  estimated  that  in  197^  the 
equivalent  of  hS?  million  pallons  of  pasoline  was  consumed  on  the 
86,000  farms  in  the  state  of  Kansas  alone  (1).  This  fuel  was 
consumed  in  tlu3  process  of  tillaije  and  harvestiat  operations, 
irrigation,  grain  drying,  marketin.':.  business  and  management 
transportation,  and  poultry  an-t  livestock  operations.  The 
production  of  fertiliser,  a large  consvimcr  of  natural  gas,  is 
not  included  in  the  above  estimates.  Therefore,  it  should  be 
clear  that  the  American  farm  is  very  dependent  on  an  available 
supply  of  oil  product  fuels. 

Common  sense  indicates  that  the  supply  of  oil  will  even- 
tually be  exhausted  at  the  rate  it  is  presently  beiitg  consumed. 
Kstimates  of  our  renainirvT  oil  reserves  range  from  30“50  years 
(2,3).  Many  ituiustrics  are  today  converting  to  electricity 
produced  by  coal  or  nuclear  powered  plants,  but  tlierc  docs  not 
appear  to  be  a methixl  available  at  the  present  time  to  allow  a 
conversion  from  petroleum  fuels  to  electricity  for  the  mobile 
high  energy  dcmar.ds  of  most  farm  equipment. 
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The  con'/eriiion  from  potroloui:i  baued  fuola  to  a more  plen- 
tiful energy  aupply  for  oven  only  a few  farm  oporationa  would 
help  the  overall  uitu.itioiu  If,  for  examiJlo,  irrigation  water 
could  bo  pumped  without  havin^j  to  deplete  the  limited  supply  of 
n^itural  gas  ( the  fuel  ured  on  the  majority  of  raidwestorn  .farms 
for  irrigation  pumps)  there  would  be  that  much  natural  gas 
available  for  other  uses.  The  simple  conversion  to  an  oleotrlc 
pump  might  iielp,  but  most  of  the  electrical  cower  gc'nerat(;d  in 
the  midsvest  still  depends  on  petroleum  fuel  supplies.  So,  an 
altorratlve  source  of  energy  that  could  easily  be  used  to  pump 
ii'i'lpatlon  water  v/ithout  deplotin.g  limited  fossil  fuels  v/ould 
bo  a better  co.lution. 

The  use  of  metiiaiuc  derived  from  animal  wastes  has  been 
proposed  as  a substitute  for  the  natural  ga.i  presently  bein,g 
used.  So  has  iiydrogon  produced  elec  truly tica lly  at  the  farm 
site  using  electrical  energy  deriveil  from  solar  sources  (and 
possibly  from  convonllal  coal  and  nuclear  fueled  gcnei-atlng  sta- 
tions duriiv:  off  peak  periodu).  This  paper  is  concerned  with  the 
evaluation  of  two  altornativ'e  energy  systems  for  pumpiiv;  irriga- 
tion water.  The  first  method  would  utilise  win.d  power  coupled 
mechanically  to  the  pump.  The  second  method  would  use  wind  j'ower 
to  generate  electricity,  electricity  to  produce  liydrogen,  and 
hydrogen  to  fuel  the  pumps. 

An  extension  of  the  second  method  could  eventually  lead  to 
jn-oduction  of  hydrogen  on  a large  scale  ai'.d  the  use  of  hydrogen 
as  the  basic  fuel  on  the  farm  of  the  future. 

There  are  mani'  approaches  that  covild  be  taken  in  devclopittg 
an  "optimal"  system  for  use  on  a farm.  There  are  ii.atiy  different 


typos  of  wir.dwvlls  currently  available  at  a variety  of  prices. 

If  hydror,on  is  f^oirr'  to  be  incl.ucled  in  the  systora,  suwo  sort  of 
stora^jo  facility  should  bo  considered  for  the  hydrO(;eu  produce  I 
dicing;  lo-.y  cnorf^jy  ders-nne!  times,  I.ittle  actual  opliwimtion 
work  has  boon  done  cn  systems  similar  to  those  consiu'-TOd  in 
this  paper,  by  cons iderin;"  the  system  as  sim.pjy  as  possible, 
lineal*  prot>ra!main.s  could  bo  ap’plied  to  nininice  the  total  cc&t 
involved  in  ostablit.-.hiu:';  such,  a system.  By  vic.vir  ; the  problem 
from  different  poi.ats  of  viev/,  dynamic  prc^jrammii-.^y , stochastic 
pro^rammin^j,  or  multi-objective  dccisioti  pvoble.a  tcctuiiquos 
could  bo  utilized.  A "syste.Ti  dynamics"  approach  is  taken  in  this 
paper.  Kathor  than  attemptii*-^  to  satisfy  coi’taln  objective 
functions  subject  to  various  cons  ti'aint.s , is  would  bo  the  case 
with  tho  techniques  already  nuir.t  ioned , system  dyramlcs  is  used 
hers  to  .model  a proposed  system  and  to  evaluate  its  porf orm.ance . 
Fro.m  tl'.cse  findin=;s,  oluanyes  in  tiio  model  of  tho  pi-oposed  rystou 
can  bo  made.  By  continuin.f;  in  this  fas-hion,  a "best"  design  can 
bo  forced  to  emerge . 


uho  syetom  uViox-  Lr.vo'jxLs^tion,  bc.'jme  n.  sii'iul’.bor  traoliig  the 
dytiarr.Lo  behavior  of  tlie  sycton.  It  v.’ac  cui.'l.'v;;  this  r:?.:r.o  period 
in  tinr.o  tiiat  f.’ortort  V/ior.or  developed  Ills  concepts  of  feedback 
aycter.s  that  wore  later  labeled  "oyberne  tics" , Sotr.e  of  the  first 
introductory  papers  in  the  field  of  feed'otek  control  systems 
(soon  to  bo  coined  "scj-voreob.anic..;;" ) wero  v.'r-ittoc  by  h-.ivo.ld 
L.  Hason.  Gordon  3,  Grewn  created  tl'.e  yorvoinoohanic.'r.s  .'.ahoi-a  b>ry 
in  the  I:  was  h>.  re  that  the  theory  of  fcvdbaok  nysto.i.s 

was  first  expanded,  tau.'rht,  recorded,  and  i ritroduced  to  the 
'.vorld  on  a lar.^e  ecale.  Dit^ital  conpiitcrs  ’.ver-o  first  used  for 
system  simulators  in  the  i?50's  at  the  Digital  Gcr.putcr  laboratory 
and  Division  6 of  the  Lincoln  Lat’oratory  uiider  tho  direction  of 
Jay  W,  Forrester,  Siitce  195b  Professor  Forrester  and  a group  of 
associates  at  t)io  Alfred  P.  Slc?.n  School  of  ;.',ana,yc:i'.cixt 

J'.ave  extended  the  concepts  cf  system  dynamics  to  cope  \.ith  tho 
great  complexity  of  social  systems  (2), 

For  many  years  tho  modclityt  of  feedback-loop  systems  has 
bccT\  known  as  "ir.dustrial  dynamics".  But  with  the  application  to 
important  areas  outside  of  the  industrial  structure,  the  name 
"system  dynamics"  has  come  into  more  general  use.  Applications 
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of  tiyutcn  iiy!\tnic'J  li.ivo  u;\dc  in  ccrporatc  policy,  uyr.niilcu 

of  '.ILaboton  n'.oJical.  synttnn,  noo.ial  I'oicoc  affeotir..' 

a b!  ic  tion. , a:'.f  tlic  boh.tvior  of  roL  O.u'oh  aini  Jnv'' lopniont 
oi-,:ani7.-i li o:-.;;  to  lint  jnct  a fov;  (4, '3,6,7). 

2 . 2 GY.iTK.VS 

Alr.ont  any can  be  T-houfjht  oi  an  a syntom,  Man  li'/os 
in  a ncrial  syntorT;.  Du-on/^h  hin  scientific  recoarch  ho  Inr. 
rncovoi-.i  some  of  r.iture'n  .'.ystenf.  r,  true  t'oro'j . 'fhronch  hin 
t’.'chi.olo^ical  brM’.\thrcuj!'.o  he  hac  creatO'.l  cov.olox  physical 
.;y  r,  toir.n  , But,  h-*  still  vlocn  rnt  fully  uniorstand  the  pr  ire  ijJl.os 
ftovorni;',.-;  tho  tohavior  of  r.yato; . 

As  used  hero,  a "cyotc;n"  is  a ^roup  of  parts  that  operate 
top.ctl’.or  with  •,  co:T,-:.on  gcal,  A oolioul.  is  a syuto  ’.  for  e.iin.’a tion. 
rancir.5  or  shivvin,::  in  a cc.y  is  a nystera  for  iCipu.lr'tir.t;  body 
tciDporaturo , Corin'.i;nica  lieu  satellites  in  spa.ao  provide  a system 
by  -which  ccmmunication  bGt\.oon  alvost  any  tv--o  points  on  earth 
is  ir.stantanioucly  I'jcsiblo. 

With  syster.c  leiiv-;  co!.':mon  aitJ  easily  defined,  '.-.'hy  is  it  that 
man  has  been  so  Ir.te  in  grnspii-vg  the  concepts  and  p'.-inciples  o.f 
systems?  Jay  W.  ferrester  b.as  poii-. ted  out  three  possible  answers 
to  this  question  (S).  First  cf  all,  u.ntil  recently  there  has 
really  been  no  need  for  unders tanJ  irg  the  nature  of  system-j.  The 
systems  that  primitive  m.an  were  exposed  to  were  n-.o-stly  those  of 
nature  and  were  .accepted  as  telr^  divinely  given  and  beyond 
comprohensior.  ar.d  control.  In  order  to  survive,  man  adapted 
himself  to  these  natural  systems  and  change  came  slowly,  usually 
through  evolution  rather  than  desire. 

Secondly,  maipv  systems  did  not  seem  to  possess  any  .meaning 


I'V  t,liuoiy.  A!',  iii.in':i  i ■)  t,  i oi!  (.•im' ■''n<t  ■lyi:  li.iiiii  ol' 

ti.uio,  i.T'oiunnic!! , uu!  I'Ollltc:;  l-i',';nri  lo  evtJlv',  ni  ui  \v.\;!  ’iriAblo 
to  il  i •;  lA  i ;li  .»!'  v>l'  the  Ini  oi.tc:  t ' oim  an'l  pr  i iii' i p 1 >;:i  diii.  (o  Llia 
tyVi'r-.jll  <;>.'m[)loxi  (.y  tlu;  :jy."  Ai’.i.i . Aii'A  I'iniHy,  i'S'ti  .ii'l.iu-  nvm 

to  r.i!ok  tlio  uri'l^'C  ly  i iw;  no  i |>li.’r:  oC  ;:y:-tonu'.,  tin  no  poini-'l- 
j)1cr.  wcro  oo  olir.oun,’  Ltint  tl»oy  wot'o  oivli;  l.oc  tod , tlurilV’,  tlio  I.r.M; 
ootUury  At  li  vo  licooim)  clo.xo  Ihit  A.ho  \i.\n-lcr  l.o  vir..'ovoi  i i;'.  ttio 
px' i no  I (i) I'll  of  rjy:’ ti':i'.!  w.x.:  not.  tho  .ib.nMic''  of  in.’r.i  1 oi'ni-i.>|ito  . 
Ifithof,  At  \v;x;;  t!io  dkl'floMtty  of  <!  .vt'ocno,  in  ■,  thnno  i c i no  i p 1 on, . 

M.'i  lii'.m.'x  t ion  invo  itily  no  l.i  1. 1 vo  1 y n,;oont.ty  ni.'';'!i'd  a lovnl  nnoix 
.xn  to  I')  .'dun.'xto  Ion  Innidlin^'-,  tiio  i.'.'oonti  i)  nonlil.)  of  •:omo 
of  iv.on'n  moo’'  oomi'l.’.'.x  n,y  ! toinn . 

l.orun'n;  f non  jixn  t I'xponi  < non  in  dil'floolt  \ ltlio\il.  a 
ntniotxino  to  Intn'ino  1 to  fno  i o -in!  ol'oi'nv.i  t Looi'.,  .loroino 
lirvinor  of  H.;i’‘v.u'd  (9)  iiao  an'und  wall  tli<i  ii.-pin' i aimc  of  ;;  tnuo  turo 
i XX  ndiioatlon, 

"dnanpi  n;':  tho  .n  t noc  torn  of  a :.nil)joi't  in 
xxnd'M'O  I xixd  t it  in  a way  lli.xt  ponnita  i;i.\ny  oi.Ixoi' 
tlilivtn  to  ill)  nolalod  to  it  no.xni  n^ifxi  i I y . To  loanix 
nAnootvifc,  in  :;ln.)nti  In  to  Ix'ann  ho.v  ttiin<’.n  ,xvo 

rolaloii f'.ood  In  .loli  i it'  tloil.  onipt',.ni  S /an.'  I.ixo 

ntniiotoro  of  .n  lobji'ot  la  pr'ol'il'ly  ovon  nonn  vatoai'lxi 
for  (ho  loan  .'ttilo  nlinU'iil  (hin  fon  tin'  .'’.il'toki  otxx', 
for  it.  in,  till'  fooit'.cr  rather'  thin  Itio  laLtnr  who  in 
moot  naaiiy  tlirown  ot'f  tlxo  trackt.." 


liYSTl’M'o  - Ol'iiN  AN!'  I'lvKDHACK 

All  nyatom.'i  ran  ho  o 1 ana  i f i nd  an.  i\i  llxi'i'  opnn  ;;ya  t>'mn  or 
ixlonod  ixy!' loi.i''. . tHo.md  !;yi;lcmn  arc  alno  c'ail.nl  " fccdh.-iok" 
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rsyrs Lciii:'. . rho  clirL'oronoo  betvv.'on  an  oprn  ;;yuli.’m  aivl  a ft-’Cilback 
syatoiii  La  that  vii  a foo(.J>i.ick  aya  tom  tk.o  b(bav\or  La  in  ('A  no  no  ml 
by  Lta  O'.vn  paat  bo'n.ivinr  'vhei-naa  an  op. -a  "yatom  ia  urnvat'O  of 
Ita  piat  tt:b:ivLor . A watch,  taken  liy  ilaalf,  ia  not  awaco  of 
lta  i iiaocuracioa  an;l  can  not  correct  itr.clf,  tima  It  may  V>o 
tb.ou;;ht  of  an  an  op^on  ayatoni.  But,  if  yon  ir.cliido  t tie  owner 
of  th'.>  watch  in  t)n'  ayatcm  ..iivl  the  o'.vner  conti  nnonaly  corroeta 
tlio  i naccurnc  Lea  in  tho  watch,  the  ayat'..n  may  bo  t;io)i;pl\t  of 
aa  a >.;lc.:od  (fcadhaok)  ayatcnn  An  0;cn  a-yatcm  and  a feedback 
ayatOii  are  nkeina  tic.rl  ly  1 1.  Ina  ta-i  ted  in  Pipparo  d.J.i. 

V.'hether  a ayatom  is  oiaessifiod  as  an  open  syoten  or  i 
foodback  system  Is  more  dt.'p'Ondent  or.  t.he  viewpoint  of  tl’.o  observer 
definjn^p  the  pin-posc  of  the  syr, tern  tlian  the  particular  assemlvly 
of  parts,  A broad  purpose  may  be  soon  av!  a Inrco  feedback 
system  consist.ifv;  of  several,  stiall  subsjc'. toir.s , oaclt  with  its 
own  purpose.  !iacii  rubsystom  may  ti'en  in  turn  be  cron.s- Idci'rd  .an 
open  or  a feoditack  system. 

Tin;  basic  structure  of  a feedback  .loop  appears  in  Fi^twro 
2.  Tho  feedback  .loop  forms  .a  closed  p;alh  st.irliivj  at  .a 

doci.;'ion  that  ooTitrols  an  action,  piassin,'^  through  the  level 
(state  or  coc.di t.ion)  of  tho  system,  i he  information  about  tho 
level  of  tho  system,  and  t.'rmina  tiarj  at  tiie  nocis ion-making'; 
p;oint,  Tlie  .available  information,  at  any  i^oint  in  l:l:ne,  is  the 
basis  for  the  cinu-ent  dec ir. ion  tl;at  controls  the  action  of  the 
loop.  It  is  thi.s  action  which  alters  Die  level  of  the  system, 

Tho  level  of  the  system  ftencrates  information  about  tho  system. 
This  information  may  bo  delayed  or  oi'''cneous.  it  is  this  infor- 
mation (apparent  level)  of  tho  system,  thoup.h,  that  is  used  as 
tho  basis  for  the  doe  is  ion  prooe.s3. 


-LT.n.'i*' .V 


(a)  Opon  System 


fi’edVack 

r " • 1 


i mni  t 


_Jj iilliFili.L 


(1))  KoodPaclc  System 


^i.‘;uro  ■-  Schomatic  of  (a)  Open  System 

(b)  FeOvdback  System 


1 
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J Iviok  loop':  ! I inlo  l.vo  ca  In,  .or.loi;  j poailivu  b u'k 

r.ivl  M'l.ativo  f'-o.lbaok , 

111  a pcialtivo  t'coil  I'a  .'k  loop,  i.h-;  t’(.’':ull  of  liio  -lolion  /'OJut- 
otoa  .Mill  '.I'O.x-l'-f  aolion,  llmo  no.  t i n,;  n p.ro','/t)i  procoaa . 
t'o'i  i l i yo  focOb'iok  loops  aro  rolilivoly  r ii'c  in  nxtiuo.  Tho 
fO'-p 'Onb  in.r  of  int.'ic-t,  <.!i  ir..'ii(.>y  in  a I'ank  'looonnt  ;;  aii  oxaiuplo 
of  pOitlivij  Xcobomk  ('a.-o  !-m  nr.'  Tho  cunonr.t  of  nionoy  in 

tho  o.coui.t  at  a ./ivon  paint  in  tiii'i.'  o ir.!  m:  no-;  tiio  'naoui'.t  of 
liU<.'."o.;  t tho  aco.aint  a.'Tn  ';  l.iio  p.font-.'O  the  aiiiount  of  itoney,  the 
coo  t tor  I ho  anmint  of  intcri?ct  oarr.ob,  Tho  ai.ioant  of  i.ioii.-'y  is 
tlai.:  l;\'ooa:c,l  ^ or;  t inuoia:  ly  'oy  t'le  i n!  or-’';; ; pai'i  aiv!  the  r.tnount 
of  ititoi'ost  paifl  inoroasoa  oont  i nnou,:.ly  ciio  to  tho  [ncrcoiin^ 
ar.o’Mit  of  r.i'nty  in  th'  r;oca.rnt.  ’.'orul  atlon  also  .’.onor  il  iy  fol- 
lows t'n ' mil  s uf  i p'oMti','  f'.'.'d  bi ..oop,  A ,:io,;.lo  coll  sr’i'.;; 
into  two  colln.  T.'iero  two  or)  Is  (h-n  oio;i  split  I'csnlt  i rt;  In  a 
tota.  1 of  four  O'''!;-.  The  r.iOro  colls  .f.o.'o  p b.  fare  'jpl  i ti.iny, 
tho  nnro  there  '-vi  1 I iv.  aft.'r  splitl.ln,;. 

In  a no.ja  L i V'.'  too  ihaoK  loop,  Mio  action  s1.riv..'s  to  raiiMain 
a specific  lovol.  A tt\o  rii'.a;:  t,*  t in  a bui’-lint'  is  a, a ox-napie  o.f 
a la-.  I i VC  fooiil'  iek  loop  ('joo  i''rpurc  .?,3.^0.  'I'li-c  .aiioeifio  level 
is  tlie  tomper.i  lure  .set  on  I'to  tiierinosla  I . If  llao  '.'oo:a  teniper- 
alurc  (i-ops  bel..'.v  tills  level,  l.iie  iv.iter  oomc::  ota  until  t no 
i'ooti  teir.iei'.a  turo  I'O.'urli.rs  the  .sp.'o  i f i oii  level  at  which  lime  tiio 
ii'calor  is  .siiut  off.  If  it  stionui  becair.e  too  wri-M  in  tiic  room, 
tlie  air  condi  tic'iier  is  svritciu'c!  cn  until  tlio  tcnpcr.atnrc  is 
rofluocd  to  tiic  sf'ocifio.1  level.  Nc;talive  feedback  loops  are 
quite  common  in  ;vt!ure.  Tlie  tlieranore/'nlu  tlon  of  tiio  luiman 
body  follows  elosnly  tli.at  of  tlie  builditVT  lliermeal.it  .already 


Figure  ?. .3.'!  - Negaiivo  KceF. back  Loop 
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<1  i r;cu2.;o<i  wll.'i  jv’rapi  t ;i  I i on  not  i:v;  to  cool  t.ho  I'Ovly  v,lion  ever 

no.J  2h' v,!r  i t's"',  nrvi  Vody  iin>  tabo  1 aotin;"  to  warm  tin'!  bod; 

wl'On  it  is  o}ii..l.'d. 

Tlio  co!)cc-j,.t  o:'  i nf orin \ t i on- I'oodb ick  l.oo]>.', , tiot.d  poo-lfivo 
and  nc-t“at.iv' , f-ov ii!.->  ',  ,-i  basir.  for  attomptinj  to  !jn  ioro  t ind  an  J 
interpret  Ihr  boli.ivior  of  .ill  ayatoma.  ’.•'oo'i'ra.'k  theory  h;ia 
rr  :co:' ' fu  liy  ’ven  o.cplioJ  to  T.i-'.oha'r  v;a  1 an.i  oicetrioal  nyntear.; 
over  t>u;  ra  il  50  yeiro.  i'ai-inr;  the  . irt  ’.'3  yea  • : t'nooe  ' mio 
princlplor.  havv?  boor,  aprll’l  to  rocial  ryatoi'.n  ij  v.oll. 

.tlOOHLS 

piode.l  la  a .in'.'r-ti  tuto  for  in  ob.iceti  i:  i t.na  t.  Lcn,  or  r-yoten 
[t  ran  le  oxprcaai.d  in  r:any  for::::>.  lionoreto  nodeli;  of  earn, 
i rpl  inor  , aixi  ahip.i  ari’  oon.rion  toyr  0.1  wi'.ieh  oliiloi'en  etton 
focun  im.i;; i n.a t ive  udvontiir'!.'-.  OoncriMe  niodo]!)  ot  varioii.'.  ob.icc  c: 
al"o  aid  n-' lontia  I a .and  on:;i  noera  in  v.'ind  tur.rnl  t-.'n  to  and  in 
V 1 17. i n^t  i5p.a.te  and  arran,.;  J^'O.ntr-  i:i  arcli  1 tec  tur  n’,  doai-tn.;;  aiid 
in  oo.ii.aanil ty  dovo’..opemi''nt.  but  oven  nioro  ooiv.xc-'.  than  conoroto 
.aodola  are  .abstract  r.iodols.  .try  cot  of  ru.loa  a. id  rol  a * lonnh  i ps 
that  doscribo  an  oiijcct,  citn.ationi  or  Gy;;to;n  is  an  abs  tr.ao  t 
mc.iel  of  that  cbjocl,  citnation,  or  systom.  Our  i.r'ntal  pro- 
I'ossos  use  oonooptc  of  a iitodol  wiiloii  wo  “vinii'ui  ito  into  now 
an  .inr.o!".cnts . Tiioao  inor.t.il  oenoepts  are  actually  an  abstract 
mcd.'l  of  the  real  system,  T.'ieso  abstract  mental  models  .iro 
.‘jubjoet  to  filtcrin.';;,  distortions,  an.i  delays  as  a result  of 
Cur  in.lividu.al  porcoptioniJ  .and  o.^perionccs.  tlioin.a  I lea  I 
simulation  models  are  a special  class  of  abstract  models. 

The  hum.an  mind  is  capable  of  absorbiivt,  adapt  in;:,  buildin,t( 
and  ucitv’:  fairly  complex  noJels  of  various  systems,  llul  llio 


I 1 p^jKinseiQiism 
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hiunan  rairJ  io  I'.ubjoct,  t.o  I'ort:'. i diM'vbacka . T'i'e  exTeot:! 

ot'  thox;'!  d f.r-vl'.Tck;!  are  ani^lif iO'i  \:hoii  the  ayatr'ii  bccorcer.  very 
cc\.ip'lcx  ^r.]  dynan.io  wi  111  f.l.vif?)  beiiavior  la  i nt.rcrtucod . 

(t)  ydn<'.al  model'.;  arc  ijfton  ill  deriried.  [ntoi'jxre  1 1 lioxi.-r 

and  aan'.ur.f' li*''t.r;  .•’bout  I'l.  il-lii'c  a i tun  tionr.  are  cotil  i nuouMl  y 
• .'hnr.'in.,';  rc'n'.ilti!!;;  in  In  v.h,.ik  the  rr.ndic  l.e  Lmjily  .ind 

orivn  roGultuv;  ’‘n  Internal  mciital  contr.adxcliona. 

( .? J It  in  ni'icit  not  nonslbio  to  reviev.'  lu'.'i'  a :nc‘:il<il  wedcl 

wae  formnl  it'id . I'-oivaceuen’. ly , aa  inniiti  ons  aro  oi'to;i  not  clo.arly 
i do  all f lort, 

(;l)  r’ar.iation  of  a mental  model  i;;  .an  individual  e.xp.n'lc.’ice 

n-.:’.' j-’c  t to  ’ n.d  1 vid'.i.il  pre .i’l ! i ae n and  peracnxl  foelirvis.  A-;  a 
ro.'.uVt  of  thin,  r.cntal.  'ludeln  are  not  o.~c-3ly  co.-.;iunicat ’d  to 
o thorn,  'f’-vo  non' ties  may  feel  they  both  have  the  same  ir.-ental 
model  of  a r.yatow  b',xt  in.  aetuality  the  t-.'/o  model:;  may  differ. 
Thoae>  xiifforenoea  m oy  bo  nlijjht  I'ut  they  car.  often  lead  to 
iieatcd  d ina.trcomonts  and  fal.ne  oonelunionn.  n 

(d)  The  c.ar.arlly  of  t'ao  human  '.T.ind  for  abnci'binp;  and 

R’ani  r'jl.atinxt  idoe.s  in  oxtromcly  v.i::  t but  it  i:;  limited.  '.Viicn 
one  l.'icn  to  r.auip'.ilato  eon'.plex  dynaxio  mental  m.cdGln,  cenfmxvon 

in  uau-'.llv  the  ronull.  i'he  ure.idcd  hu.aan  rnitid  in  jra  t not 
.adcijuaro  for  conn  truotin.^  and  ma  n i :)u  I a t ■ n,t  ifyn.’mic  .aiodc),:’. 

f’yntcm  uynamicn  provide  a foundation  for  lrar.nl. a tin*  montal 
modelr,  into  i.ia  cho'x.a  tical  modol;;  v.-hich  o.ui,  wi  tli  the  aid  of  a 
computer,  be  manipulated  with  relative  oa:;o.  by  tran'n l:\ti 
mental  model::  intx;  .m.atliem.t tic.il  moxbrl:;  mofit  of  the  proble.ms 
.an.acoiated  with  montal  inedol:;  are  a.lleviatcd.  With  all  modola 
exprensed  in  the  :'.ami:  tcrminolo.yy  vienfnnion  about  term:*., 


.arrjump- 
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,\onr,,  ar.’!  oonrunic.it Ion  of  motlcL^  i'J  rcfUscod. 

.’loiluj  : i'linuld  not  b^.*  judged  a,';’nLnot  ;r.ia.~in.iry  porf'-'c  tlon. 
.father,  tii''y  shculc;  bo  juJjod  in  tonns  of  lio.v  well  tlioy  describe? 
the  Gyat.'i'i  coir.pired  to  exlntln^;  ".sodelc.  The  cor  tain,  ty  '.vith  (.■Iiic'n 
Kodoln  cho'.v  t)ia  cor'-cct  time-vary ir^;  results  of  ^hc  noeb'!].  stato- 
ir.ents  conparel  to  t’.io  same  c.r'.eluoion  reached  in  oxteudin,'! 
inontal  models  i r.notiior  base  on  v.t'.icl?  to  judge  Ihv  !r.at!u?ria  tical 
ir.oJcl.  y.oJels  Gh.ould  tf’ien  bo  judged  not  on  an.  abGolute  scale 
that  upu'auds  cr  co.nde.ru’.G  them  but  rather  on  a relative  scale 
that  approves  them  if  t)icy  sa  f.icractorly  clarify  the  kr.own 
facts  about  a sj'stcm.  One  model  is  better  than  anotlior  only  if 
LI  better  cvorc.unica  ten  tl'.o  mo^'.tal  and  pliyeicr’.l  s ta  tenents  of 
t!?e  .system  with  a l\i.gh.er  dogreo  of  reliability. 

2.5  SiXULATlOH  TiiCii.ll'.lUaS  v.5  .ANALYTIO.-J.  SOL'JTiONS 

Simulation  is  an  (?/perir.'.onte  1 attempt  to  represent  tlio 
behavior  of  system.s.  ft  consists  of  a s ti  p-by -,:tep  apf  roach 
to  "simulate"  the  dynamic  (time  varylr.,;)  boh.avi  or  of  the  real 
syr, ten  th'  model  reprosontG.  An  an..lytic  •'olutlon  expronaes 
t.hc  system’s  conditii>n  in  tci'ma  of  any  futut'e  time,  not  just 

in  tenas  of  short  time  intervals  between  comrutations  .a.5  in  a 
simulation  approach.  In  addition,  the  form  nf  an  analytic 

solution  can  tell  much  about  the  general  n.rturc  of  the  rcspon.jo 
of  the  system  without  havin,';;  to  carry  out  aiii'  numei'ic  com.putatlons. 
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To.'  <!  ■.  iroplc , conriirlci'  n 'j iwpi.i I'lod  invontory  contcol  :;yjLoia 
( n^oUMir'.'  I'.o  dolay  bciv/ocn  order lii,';  of  ^^oocio  a.oi  tiv.-ir  receipt 
Into  I nv  Cl'.  ;ory ) . '■'i.'^.uro  d.5-l  shews  tiic  siitiplo  f it t -order 

r.Cfptlve  f'cdbtc.k  loop  tliit  r.cdols  the  systoni.  The  of 

the  r.ystcrp.  modeled  is  to  .iiiir.tala  tl'.o  desired  inventory  D[ 

■.vhioh  is  a coivstrtrt.  I'rsu.ne  t.h;’ t i..  re:  liaridi.se  oan  b->  eiUior 
ordtre.i  or  returrwi  to  the  r.ariuTnc lurcr;  that  is,  the  order  tale 
OR  car.  be  eitli'.r  pcsilive  or  nei^ativu,  fa  a t tei.io ti nj  to  brins 
the  aein.al  inventory  to.vrrd  th.r  dosire  i inventory . ti:e  order 
rate  itiist  inovoasc  positively  as  i r-.visnt ory  falls  a.'.d  noftalivoly 
as  invcr.tory  ireT''a.'!os  beyor-.d  the  desired  inventory. 

Assuming  a i Ir.onr  relationship  for  the  order ir.g  rat-.j,  the 
syst-’.ii  could  be  deoerib.id  ir:  tiiuntion  fo^.i  asi 

OR  - (01  r) 

who'.'c  OR  " ofdcr  rate  (uni t.a /week) 

.AT  - ad just.i'.i'iit  tiiro  (v.-eeks) 

Dl  - de.sirci  inventory  (units) 

I -•  i Inventory  t units) 

Tilt?  ad  j'u3  tment  tine  .AT  factor  is  included  to  .make  the 
equal  i.en  il itre.nsi on  t J ly  bala.nced  and  by  being  a constant  cst- 
allirdu’s  the  lii'.earity  of  the  system.  Assume  that  the  de-sirod 
inventory  DI  is  6C00  units  arvd  th.at  the  adjustment  time  (the 
time  that  ar.y  cin'i'ent  order  rate  would  require  to  correct  the 
inventory)  i.s  5 weeks.  The  equation  describing  the  system 
then  boconcoi 

OR  ^ i (6000-1 ) 


I 

1 


Figure  2.5.1  - Simple  Invcntoz-y  System 
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Assiiiiiri,';  .-.liio  that  the  ir.Uial  invoatory  io  tOOO  unlta,  tho 

iniilal  order  rate  hecoinoa  1000,  How,  if  thin  order  rate  la  in 

effect  for  ?.  v.-eeko  'bafora  i new  order  rate  ia  cal'  iila ted , th.cru 

j will,  b?  2000  ur.it.3  added  to  inventory  which  tiie.n  iecomoa  jOCO 

! units.  Thia  now  value  of  I results  in  a new  order  rate  ot  600 

i 

units  per  week,  by  contiuuifv?  in  this  fashion  it  lr>  possible 

1 

; to  ci.nulate  the  dyr.a.Tic  be!\nvlor  of  the  systo’.n  through  t.’.'Tie. 

I 

Proceeding  v.ith  tliese  caleulations  yields  the  results  s/iown 
' in  Table  2.5.1.  "'l-.sce  rcsslts  era  slu/.vii  plotted  In  Figure  2,5.P. 

i This  demonstrates  tho  siaiulation  approach  to  a solution.  Mew 

i consider  the  analytical  approach  to  the  sam.e  pioblem. 

Duo  to  the  relatively  simple  structure  of  this  system, 

! 

it  is  possible  to  arx'ivo  at  an  area.lytic  solution. 

■;  CR  - an-  (01  - I,,)  EQ(2-l) 

j n .r.i.  fl 

^n+l  " Interval)  EQ(2-2) 

I ■ ‘ 

; wiiore  the  subscrips  n and  n+1  de.aote  the  present  ar.d  next 

i 

t.ine  interval  values  respectively.  /Issunir.g  that  the  rime 
interval  between  successive  siraulacion  computations  has  become 
■ v.ai' I shingly  smallj 

twl  - In  " O'  " ri'  <='  - I„>  01  '00(2-3) 

K “ -A'  (”'  ■■  O 10(2-4) 

lit  o Af  «>  “ aV  (»0  I'«<2-3) 

; Equation  2-5  is  a first  order  differential  equation  whose 

solution  yieldsi 
f 

I 
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r(t)  A I- 


KQ(2-6) 


As; r,u;ni r.''  uri  ir.vciitory  of  it  i'.;  pDiiolbl.o  1.0  <.  jtco.'ii no 

the  value  of  oon:;i  .u\t  A^, 


Thorafofo , 


1(0)  ^ A^  + r)[ 


. [,^  - ni 


I(t)  (I^  - + 


KQ(2-7) 

Ea(2-.'3} 


Subc  ti  tut  Ins  in  tho  ccnstan'c  valutu'  annuin.ioJ  in  tlio  s;!.Tiulat.i.cn 


loiution: 


I = 6000  - 5000  0 

By  su'crrlitucins  any  r.lno  (in  •.vookn)  for  t,  t.h  ; invoncory  at 

timt  tir.o  0:in  bo  00;  rou  tod  ciirocMy  wif.hout  c;oir..';  thron.jh  tho 
ntep-by  -n  top  calc-.iia  tiunn  requ.w'cd  in  ihc  aiv.ula  tion  uvonoduro . 

The  negative  exponential  nature  of  the  ayrite::!  shewn  in  the  plot 

of  tho  nianlaticn  solution  (P'Lgure  2.5.2)  is  obvious  fiO;7i  the 

form  of  the  analytical  solution. 

Since  tho  analytic  solution  of  a synstonys  be)iavior  contains 

so  much  infori;'atio;i  arid  loinco  it  allow::  dinect  comput.ition  of 

tho  condition  of  a syGto:n  at  ani'  ;?pec'Lfied  time,  it  :nig!it  bo 

prcGum.ncd  that  an  a;:::  lytic  solution  should  iilways  bo  obuainod 

for  every  system  under  study,  however,  this  is  not  pos:>iblo. 

rrocent  mathematical  urccodures  do  not  permit  f oinnu l.a tVon  of 

an  analytical  solution  for  many  complex  systetas.  When  an 

analytical  solution  is  beyond  the  scope  of  today’s  mathonatico, 

simulation  techniques  are  used  to  model  th.e  dynamic  tcliavior  of 
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the  system. 

‘I’liO  u;;t‘  or  r, i ,-m l;i L Lon  heon  in  th(i  i«a::t  hue 

to  the  cootfj  on!  tline  inviilved,  Ilavifj:'. tion  table. ; nilin.';;  l’r.'0..i 
the  toi'O':;  tli'vv.  weno  coinputou’  untti;;  .n  inu  lat  ion  te.^hniiiuo;; 
rhow  the  rolati’.ely  Ion''  hinl.ify  of  the  iiru.l.n  tion  te'';hni'iuo 
(3),  It  not  ur'.til  alK'r.t  I.9j5t  \ih<:n  co'.ipu  tati  on  rU-nt 

be  carried  out  on  1 co.-nputero,  tifit  the  tine  .'’nJ  corjtn 

n.'jcoc  Ir.  ted  with  -ei'i'.ulatlen  teoluilcuen  '•■'•■■vo  reciucoii  to  the  point 
w!iore  a i "rj la t j.or4  nolnttcn.:  wore  not  o.'ily  mar.’  ca''.ny  roaolicd  tlur'i 
.antlylic  aoluLLo.n..  i'lit  were  i ne.xpen.a ive  .rs  well.  Today,  with 
i.t\o  aid  of  hip;.h  ."reed  eleeti'onie  oonpu  lera , a len.^LI'.y  sirn.il.ati  on 
of  a complex  nystet?.  t-at-cr;  ot\ly  a ftn-/  nooonds  anil  coo  to  only  a 
few  cents. 

3.6  I'LOV;  ir.\GR,i:.',3,  COtrUTATIONS,  .MJD  !''.rJ,\TI0N3 

i’ \ "'I  picture  is  -worth  a t.hounanil  wo'-Js”,  or 
equal  ions  as  llie  case  may  be.  Often  it  ,i.s  easier  to  jjrasp 
the  idea  of  the  v/orklnt's  of  a system  if  t!iat  system  can  be 
disirramed  effectively . Thac  is  the  ifoal  of  flow  dia.qranii n;;. 
Whereas  the  equations  of  :i  system  focus  on  the  composition  of 
the  elements  that  make  up  the  system,  the  flow  dia.qr.-’.m  rhoul'l 
.slio'w  how  these  oli-'monts  are  interrelated  to  include  feedback 


loops  and  how  they  effect  ttie  system. 

Tiio  industrial  Dynamics  Kcseareli  Group  at  tlte  M.T.T. 

Sloan  School  of  Mana.iemcnt  lias  developed  a set  of  .slandaid  flow 

tiiaijram  symbols  for  dynamic  models  (10). 

l^aVKhS  - Levels  are  represented  by  a roe  taiytlo . Calculation  of 
all  levels  involves  integration  of  incomlnq  and  on  hqoitvi  data 
that  are  controlled  by  rates.  The  letter  catnip  of  the  reja-esonted 


poi’Uon  or  the  sys b.'i.i,  tho  full  namo  of  t;.?  ropoo ■ ciUcd  portion 
of  tho  oyjto.’i,  nnJ  Iho  rqu'ttion  nunl/or  ao  a crooa  rofororoo  to 
tho  forwal  'Icrinltion  of  tlu:  cyat'i.i  I'.odcl  aro  l.Kl.udod  in  tlio 
roctary;lo  ( aoo  ? 'A  .6.1)  • 

r\ATI-:S  -•  Katu  p.iati.o;v5  of  a oyatoM  aodc-l  aro  tho  policy  rtxtc- 
uoru;;  and  doflr.o  the  t'lc>v  atz-o;.,-;;  i fz  the  syatem.  'ilio  only  inixit 
to  a rate  la  iat'oi'rr’.atl on  a.rid  the  n;!ly  output  la  in  tho  for^)  ot 


coa troll  in..;  riov.-o,  A ayir,i;o.Uc  v.i.lvo  la  uacd  to  rer.r'oont  a 
rate  due  to  th’ia  type  of  "controllin.y”  output.  The  lottor  group 


of  the  reproaonted  portion  of  tlie  ayat.cm,  tho  full  mine  of  the 
repreaentod  portion  of  t'r.e  ayatoin,  the  equation  nurabor  na  a 
croaa  re foreiV-'O  to  the  formal  nyatoin  model,  a.vi  the  infermition 
.inputs  cn  -.vri.ici'.  the  rate  uopenda  should  h'>  aho'.vn  with  the  rate 
symhol  (sec  Pigizrc  2.6.:?). 

auxiliary  VA.l.IA.',Li:s  - Auxiliary  variabloa  lie  In  the  fh-.;  channels 
between  love.'. 5 ind  I'.atea.  They  are  actually  pnra.s  of  the  rates 
but  are  subdivided  separately  bocaviso  they  represent  ir.d- p'eo.uently 
moaning  coiicepts  of  interest  or  because  they  maite  tl'.o  comp- 
utrrtien  o.f’  rates  Ics.?  complex.  A circle  is  the  sym!)ol  I'oc  a.n 
auxiliary  vari.iblc  and  should  izioludo  the  ablu'cviation  of  (.he 
variable  n.im'',  I.Jio  v.criable  name,  the  orpiaLion.  nunber  as  .i  cross 
reference  to  the  fci'ir.al  model  of  the  system,  and  the  inpjt  .and 
output  flows  ( SCO  I'vgure  2.6,3). 


I'l.O'A'  J.INbS  - Six  typ'GG  of  flo'.vs  are  shewn  in  Figure  2.6.6. 

Use  of  the  appropri.ite  flow  wilt  climirate  the  need  to  l.tbol 

ind.ividunl  flows  (with  I'e.spcct  to  type)  \/itlun  llio  flow  diagr.a.a 
of  tho  system. 


INFOilFATION  TAKF-OFF  - A small  cirolo  at  t!ie  information  source 


(D0)3- 


budf'ot 


2.6.3  ■•  r'lov/  Oij^srram  for  aTi  Au>  i. liary 


[nt'ornia  iloii 

Ma  tcz'ial 

Orders 

iVcncy 

Tooplo 

.■-'quipii'.enc 


* ...  - - iv 

.S 

'1.-, 

>s 



Fizrure  2.6.4  - F\.o\i  Linos 


and  the  InfhfMation  flow  line  repronanto  information  tako-off. 

This  '.'epv'OGen.to  the  re.eoval  of  Ijiformatiou  about  an  oiornont  of  the 
ayctem  '.’.'ithout  affectin’;  that  ‘.’lomen'Ci  In  goro  InGtanouG,  though, 
tl>c  information  taken  froiii  a flov/,  rate,  level,  oi’  a’uaL.lLary 
v.'ill  he  UoC'l  to  charge  ’tiiat  flo".’,  x'ato,  level,  or  auxiliary  (nee 
Figure  2.6.5). 


dfSIcM  PAFAMETFRS  AiS  GOhSTAhTS 
conGtantG  do  not  ohaiv^e  within  a 
tet'.veen  succesGivc  a inu lati ona , 


Syotor,’  para.motorG  and  other 
given  Ginulaxion  but  may  charge 
'Figure  2.6.6  Ghewa  the  ayrabol 


for  a constant. 


SOliRC-FS  A.'O  SINKS  - In  Gome  cases,  the  source  or  termi 'nation  of 


'lo'-,'  has 

:io  ir 

’ifluensve 

cn  the  sysl; 

m itself.  In  tills  case  it 

assumed 

that 

•uho  oO’': 

rco  or  sink  i 

.s  "infinite”.  Sues  a sou;n 

never  be  exhausted  and  such  a sink  can  never  become  full. 
Figure  2.6,7  shows  tiie  symbol  for  a source  and  a sirJ-;.  .'luch 
symbols  have  no  dynamic  characteristics. 


OTHER  SYMEOLS 
Oliver  sr;ecial 


- Logic  functions,  delays,  table  functions,  .-jnd 
functions  hvive  special  flow  diagram  symbols. 


These  are  explained  in  further  detail  i’n  Appendix  2 as  well  as 
in  (8,10,11,12). 


Cornutatlcns . Simulation  is  a step-by-step  co-mputa tion 
based  on  a sot  of  equations  describing  a given  system.  This 

s tep-by-step  procedure  calls  for  some  sort  of  dcfiiiite  computing 
sequence.  System  dynamics  models  con’tain  oi'.ly  two  basic  types 
of  equations;  levels  and  I'atos.  Other  typos  of  equations  include 
auxilary  and  supplimc.itary  equations,  constants,  and  initial  value 


coi;-.vr.\riT 


OX' 


CON^jI'ANT 


l.'^nre  2.6.6  - Flcn;  Oi'igr;i:n  Tvir  a System  rai'ui.ia  ter 
or  a Sons  t int 


Fi?;nro  2.6.7  - i-'.low  Dixirram  for  i\  Source  and  a Sink 
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ir;ual.i.oru:.  Tho  tirno  sii'.ii  DT  u'-.ci-l  in  .■.vryitv;  out  t!\o  ntcp- 

iiy-'Jic]}  ca  ionla  ti cr.r:  of  a ‘Jirnul  itiori  i.;  ilno  vory  i 'iip-ji' Lniit , 

‘i'lrac  nan  bo  br  ■ ron  into  tiiroo  rl  in  tin.c  t nor.  torn ; pint, 
t,i-eccnt,  and  futuro.  Tlu  timo  iioiatLon  ccnvintion  utiliv.od  in 
nyslon  c!yna.aicr>  is  i .1  ins  ti  aicd  in  Fif^ure  Z,6,'},  'i’iio  "pioscnt'* 


instant  in  tlwc  is  d.  nip^nated  by  'K',  an  ir.stauL  c.-o  ui’  time 
intorvat  in  tl'.o  past  by  'J',  and  an  instant  one  TV  tlno  interval 

in.  tliC  I'utura  by  'L'.  Values  I'jr  levels  at  'K'  can  be  <;  ilculated 
from  tiicir  values  at  'J'  (i'rorr.  pfei'ious  ca ' eui  a t ions ' and.  from 
the  valvic  of  fne  r ates  involved  darir^  tiro  intcovai  'JK'.  Given 
this  infor-mali on,  the  values  of  tlio  rratos  durin,!  tin-’  next  tine 
period  'KL'  can  le  (ra'ieulated,  Tiiia  in  the  situatioa  slsown  in 
Figure  2.6.d, 


'I'hc  sequel: 

.ce  of  calculation j 

•cay  now  be 

vepocted 

by  advatic  1 ^’^T 

the  timo  indii'a 

•torn,  J,  K,  and  h, 

by  0 no  01’ . 

This  is 

shown  in 

Figure  2.6.9  . 

All  c:alculatio:::; 

arc  carried 

cut  exac 

tly  the 

same  as  bo  fore.  The  ti'.r.o  interval  that  v/as  pr-cvicusly  'KL'  is 


no'.v  'JK'.  This  prooesn  is  continued  until  t'ne  sirraili tion  has 
been  carried  as  far  an,  desired. 


Fjjk! tj_ojrn_.  The  variables  and  cenraunts  of  a system  are 
i'epresGnted  in  equation  form  by  sy’i\!bols  (or  aiibrovial  ions ) , A 
standardinr J style  of  cstabllchir.c  those  abbreviations  is 
essential  for  practical  purposes.  A symbol  (or  abbreviation) 
for  a varl.ible  cr  constant  in  a r.ysteir  consists  of  up  to  seven 
characters,  tho  first  bciiu;  alphabetic.  All  variables  are 
followed  by  a period  and  a time  postscript.  1 cvols  and  auxiliary 
variables  carry  a sin.tle  letter  postscript  indlcatin,>;  the 


h'P.JK 


1.2.  K 


r*"  r\ 


R2.;\L 


.'’.6.3  - Computation  Soqucnco 


Flnui-o  2.6.9 


Computation  Soquproc  I'l.'ctondod 
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point  in  tlxo  for 
letter  pootocrlpt 
apply.  Cor.r.  lant;; 


v.hlc'a  t!iO  value  appllon,  ilaten  carry  a two 
inljcatin,^  t>u:  tine  interval  over  which  they 
are  clwirac  toriaC'l.  by  no  poc  toci-ipt. 


an'/lifi  bQuATICii'l  -•  A lo'/e  1 equation  leprowent'!  a re’uorvoir  (or 
acoraulator)  tbnt  -ocurula tea  cho  ratea  of  flow  that  flcyw  Into 

and  Out  of  th.o  reservoir.  \ nev/  vaiuc  La  calculated  by  nddlrv; 
to  or  ctibtr neti:t;  from  th.e  previcua  value  the  iMowc  chat  have 

occurred  durln.q  the  D?  tli.'.e  interval.  Ihe  level  equation  format 

is  given  below. 


.•hero 


L . . j i-  ( D? ) ( ;i  I . jk  -rout  . j;<) 

L - levol  equation 

L.K  - value  of  L corepu  ced  a.t  K (units) 

L.J  - value  of  L fr'^m  -previous  time  period  J (vi'.lts) 
DT  - length  of  tine  Interval 
RIN.JK  - rate  into  L during  ir-itcrvol  JK  { ur.i  t.e/tino ) 
ROUT.JK  - rate  out  of  L du'cir.g  interval  JK  (units/tino) 


The  level  equation  is  th.e  only  equation  type  t.hat  contair.s  the 
time  interval  DT  notation.  Arov  number  of  rates  may  be  added  or 
subtracted  fror.i  t.he  levol  at  one  time. 

The  levol  equaticn  performa  the  process  of  integration . The 
above  equation  could  be  vritten  in  the  notation  of  calculus  and 
differential  equations  as 

b = Lq  5o 

Level  equations  are  represented  by  the  letter  "L". 


r 


H/'v'L’E  E';UA':  ; .'"iS  - Kat:  dafia.o  !io  t'lo'.v;'.  within  a nyote'.i 

ax'c  c^ntnol'.od . hy;;  iovc.l:;  a;.'!  con'-tant::  ju'cvid;  ciu  input 
tc>  rnto::  i’v.  r ito;;  in  ti’ni  contioi  thu  fLC''.i  L , t'ron,  or  t>;  tv.-oi.'ii 
iovoln.  iiio.  f orn-.c.  t op  th,‘  .-ato  lion  in; 

!i  a.KL  i(]t>vols  aivt  c.  r.nt  xr.ln  o'"  t]v  nynt’nii) 

Rato  c\uavict'.3  aru  dorotoi  ty  ;h.;  loiter  ".V. 

A'JX’I  nr.ar'tr  I-.QL’A  I’.'cr:.';  - CPton  mto  e'ln.iti cr.  will  be  '-.liX  ividod 
to  enh.r.  nc'.'  ito  r.'ai'.in;  or  to  oapr:i'ato  aa  i:itor:.:od  la  to  value  of 
ir.tcrorft.  Thin  in  aoaoxpl i nhed.  by  the  uoo  cf  auxiliary  o p.iatiortu . 
Sinoo  auriliaxy  cqu'.'riors  oxrr  .'■.r;  tiie  value  of  a prrticular 
qu'.’.tity  at  a poir'.t  in  tiT.:,  ti'.oy  are  followed  by  either  "J"  or 
“h”  tine  sub,:oript3.  ’■'.he  ure  of  euai  i inry  equation:;  in  illustrated 
be  lev;; 

R 


CR.KL  - ;^(ni-i.x; 


A I , . K/CC 

A Q.K  EAl'/RP.K 

Auxiliary  e.puitior"  are  do;; icr.ated  by  th.e  letter  "A". 

fil'Pi'l.iy.EM’ARY  " ti'JAI'IPRS  - Sup'plimentary  eu  ua  tiono  define 
v:'.riableo  ti-.at  are  not  actually  variablLU  in  the  aye  ten,  but 
v.'iiich  coTitain  inferr, avion  wh.ich  it  is  der>ircJ  to  have  printed 
cr  plcLtcd.  Tho.oo  variablco  are  the?)  urrd  only  in  prii'itijv? 
and  plottin.j  ir.ntructions . 


IMTIAb  VAI.UE  Eqt’A'i’IOry  - All  lovoin  must  hive  tiioir  initial 
values  specified  before  a sinul’tion  can  be^in.  Th.oso  values 
arc  ncoocsnry  for  do  termini  n.q  the  flow  rates  over  the  first 
time  intci'v.al  (0  - DT).  These  equation:'  are  designated  by  the 
letter  "N”. 
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Till’  iiiMw'-rU-  i:;  U\i->  li.uik'’ t.  oT  tin'  nriU"'!  ;'ii  ili';i,  iivl 

(.Ik!  ilnl  toil  i o,  llio  liroii'l  I'.rik.,' t of  i lie  wooM.  'I’lu;  i.onloi'ii 

jilalioi  Ilf  fill’  iiiii'wool,  ;\-(i  alao  o'.oh  in  av.'i'  ' iMi.'  wiiA  Ci'.oi'cy,  It 
.1:;  In  t’li:!  afra,  Hu-n,  flial  Uio  I’li'pv'aoii  wint  I'owovoi!  i o:' i o.a  ta  on 
riyntoin'.  (V/fl  ) ua'  oi’naiilofcii , 

v.'ltid  ilat  i (T'Oi'i  tlio  Woallioo  liui'oaii  i.'l'fi.oo  in 

OiHlfio  City,  K.iti'ia:;, 


i :!  u.aai  in  I ho  auiiulation  of  Iho  i,',a’  I'roia.'iio.l 


Uiiii.i . wind  d.it.a  mnnl  t ocuuij-!  L t-  tor.Jn  Lo  bo  a " i.y  pir;i  i” 

year.  'L';i  ; nei'.ir.er  i’oaUi::  (luri(',_;  irrif.atvon  wi  ll  'o-> 

I'laoo  h.-...:  u'/era  ;(!  '.'Ir'.i  ei'iocbitlos  vnlliin  onc>  lialf  a a tat.d.u.  i 
• i of  a <'1')  yc'u'  ■wceo,^':  Ccc  tlio  nrea.  'I'ho  eLhev  nuw 

eenh  have  avoiu.p'!  wia.l  ve loe  I tiea  n Iv/o  :: ’.anJard  ' 

(levial)  oiiJ  of  a •’’di  yea"  avorn.a;.  '.Vic;.!  d-.  ta  v/ero  reaoj'd.j'i  every 
3 hour*:.  '1'!'..:  data  v'a;:  ohta.i.iied  aL  a oi  u ii;elcrn  ahovo  the 

I -.re a’  n. 


'Ihe  WPI  syr. 

toms  a 00 

d 0 s i 

L.pied  to 

provide  ei'.or>  y to 

IrriciLe 

160  acres  'with  6 

to  8 ao 

ro  -Ir 

'.dies  of 

water  per  montit  di 

i ' 1 11.3  June 

July,  and  Aeyust 

, Tt'iC  V.' 

at  01- 

t*'ble  is 

assuiii'.d  to  I'o  at 

2 CO  feet. 

3.2  mop:;!,  ^TducTUhi;  - ottcon  i 

!'■' '.13 liSil*  :-,v.  Lie  eT  v''pt5.oii  I ap-fea’.':'!  i n Pl,;\!VO  3,2.1. 

In  I'.liL..!  i'.yalcM'i,  Ihe  .■.'.in':  '..ueu'.;  a '-.'ind  tui'blne  (wi  ndiii i .1 L) 
v.l'ic'h  io  coupled  lo  u rly.diue).  throci;;i  a /'cav  box.  The  onoray 
jii  tc.e  ic  cilrecLly  ce1.;iitd  to  vho  velocity  OL'  the  ■'Ir.d. 

P'ut  the  .vlrvi  dooc  not  blow  at  a steady  cprod  for  more  th.''n  a .few 
f.ccondc  at  a tine,  do  !ho  flywheel  ac  lO  ao  a buifoi'  to  -.ven  out 
the  hi^'jh  and  le\.'  wlr.d  iMfio..!;;  and  It  nlco  neevee  to  n tore  a 
oortvaln  anour.t  of  on":i-,,y  for  uno  when  the  wind  i ::  not  blowln^. 

'fho  onore;y  in  the  fly v;lw;ol  l-j  icnred  throe, 3)1  a traravilcoloii 
to  r mechanical  pr,’..p  w'.iero  I’t  door  liio  work  of  pu..iplrv]:  the 
lrr5,n'itlcu  water.  A farn'.  trae  tor  provides  a bn.cJc  up  for  tliooe 
periods  wh.cn  the  win.!  la  not  tio-vln;,  and  tho  f.lyi.l.eel*  0 tiu'oil 
ener,';y  In  low.  I'ach  component  In  the  r.yntoM  has  its  own  ef- 
iMoif'iicy, 

A casual  loop  duu’;vam  of  Option  I is  shown  in  Ki,;uro  3.2,2 
and  tho  DYMIMO  f.lov;  diagram  of  the  mode.l  is  rdiown  in  I'i^uro  3«2.3. 
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Option  I 


1*  If;! I i*f.' 


^wind  speed 

•Sr 


. ■*  f'otor  size 


capital 

cost 

($) 


RWS  - rated  wind 
peed 


CIS  - cut-in 
*7^  speed 


wind  speed 
ad  justed 
( knots ) 


rotor  capl  tal 
^ cost 


* ^ ■ height  > 

t. — scaling 
ESWZ^T  'k  I'actor 


'power  fro« 
the  rotor 
, (kW) 


' ener^ 
surplus 
wind 

i2nd  index 


RE  - rotor 

efficiency . 


power  to 
(kW) 


total  energy 
surplus  wind 
2nd  index 
(kWh) 


F'WEC  - flywheel  / 

■TT  energy  capacity  / 


the  flywheels 


total  available 
energy 


Fi^^FP  - flywheel 

efficiency 

' iL2^“  Gear  box 

^ ~ efficiency 


' energy 
surplus 
wind 
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Figure  3,2,3  - Flow  Diagram  of  Opt 


r.jvol'o.  'Jher'i  are  nix  Xovola  in  Option  1.  'i'oial 
o.-/nLl.aX'Ie  encvr/j  i'Ai  in  trio  total  amourit  of  t,'r:or(;y  av.iiiablo 
an  any  point  in  tine. 

Total  cnai-^r/  a’Ji’plus  ( .'i.'vJ)  lot  i/idox  TESV/l  is  a level 
which  accii.Tiu la  to ■;  aiid  kcepa  tr  ick  of  wlivJ  enerjry  tha  t is  avail- 
able from  the  rotor  but  can  not  be  put  to  uco  i;r;necl lately  ox' 

.'j^aro.l  in  the  fly./hec-l  fen-  future  use.  This  surplus  energy 
cccui'S  when  the  flyuhoel  scoxar^e  5s  fall  ljut  the  is  still 

Mov.irg.  TSSV.'l  is  used  in  e/ali-atiry;  tl'.e  rotor  anl  fly./heal 
coDibi.'A Lion.  A large  T'lS'.ii  v.oi  J.d  indicate  that  the  rotor  is  too 
large  or  that  the  fly.vhGol  is  too  small. 

Total  energy  surplus  (wind)  knd  index  TESW2  is  a level  which 
nccur.ulatos  and  keeps  traeJe  of  v/ir..i  energy  that  can  not  bo 
gc.i-.ratod  duo  to  the  Tate-i  wii.d  speed  pv.'o  of  the  ays  tsu.  liuppcso, 
for  exan.pl G , ciiat  tiie  rated  wind  speed  of  the  sysxsiii  is  lO  knots. 
Th.c  energy  avaiiahlo  from  the  v.'ir.d  is  directly  pi'oporrionnl  to  the 
cubs  of  the  wind  speed.  So  as  the  wind  speed  ir.creasos  fro-n  5 
knots  to  7 Irnots  to  10  knots  xhe  system  will  bo  able  to  extract 
more  and  more  energy  from  the  'vind.  But  for  any  wind  speed  greater 
than  10  knots  the  energy  extracted  by  the  sys  ;cm  •.vill  be  the  same 
as  at  1C  knots.  Thu;--,  v/ith  the  wind  blowing  at  12  knots,  only  10 
knots  v/ortii  of  enoi-gy  can  do  extracted.  TiiS'A’k  keeps  track  of  the 
difference  between  this  amount  of  energy  and  the  amoxint  of  cnex-gy 
available  were  the  .system  able  to  extract  all  12  knots  'worth  of 
energy.  If  simulation  of  a system  results  in  a l.xrge  TkSWd  then 
the  rated  '..rind  speed  KWS  should  be  increased.  There  is  a limit 
to  how  largo  RV.'S  may  b>":eme  dopendiiup  on  the  system.  Suppose  the 
system  has  a mechanical  device  coupled  dix'cctly  to  the  mrchanic.il 
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outii-nit  cT  r.ht;  i-otor  at'.i  tliat  tho  rotor  io  rat(;il  at  10  knotn. 

Aijku.t'O  a.).:o  that  the  rotor  bo^jin  to  c.;tract  o.oor;"/  t'ro.n  the 

wlr.a  at  .'.'.out  4 knots,  k’hen,  at  4.o4  knots  iho  mechiinical  Ot-vice 
coaplcd  to  the  rotor  .-/ill  be  opera  tin;j  at  only  10^  of  its  ratir,^. 

If  RVt’3  3.S  increased  to  R1.54  knots  Lhen  the  dovice  v.'ould  bo 
operatin/i  at  only  1;'  of  its  latiiv;  at  4,64  )tnots.  This  is  v/liy 
RW'j  must  be  United.  A ranp;e  of  10  to  15  knots  between  the 
speed  at  which  the  rotor  begin;;  to  exti'act  oi;er<:y  and  the  rated 
wind  speed  is  approximately  as  large  a;;  most  ir.eclianical  devices 
will  tolerate  ar.d  still  operate. 

The  system  described  by  Option  1 could  theoveticriliy  pump 
irrigation  water  24  hours  a day,  every  day,  for  the  xhree  sumner 
laonths;  Jur.o,  July,  end  August.  2nt  such  a uonstraino  v/ouJ.d  not 
correctly  describe  the  pli^'sic;-?.  sltuution  ar.d  would  coat  a 
considorabl>;  amount.  If  the  system  ;vere  to  be  pu;.'?pin,g  24  hours 
ever/dr.y I more  v,-ator  "/ouJ.d  to  av.eilr.’.'lo  than  ’would  be  n.eedfed. 

More  realistically,  the  pump  v/ould  be  shut  off  durir,g  tiir.e.s  ’.vhen 
t];e  ’//ind  were  not  blc.v.irtg.  Those  tines  'when  li;e  pump  is  shut 
off  are  referred  to  as  "do'wn"  time.  And  the  daily  amounts 
of  down  time  arc  stored  in  tlie  DOV/N  accumulator.  The  amount  of 
down  time  allowed  per  day  is  upacified  in  the  program  as  down 
time  allowed  DTA,  nhTOT  accumulates  ti;a  total  amount  of  down 
time  occurring  in  t)ie  complete  simulation  run. 

When  the  system  has  already  been  "down"  for  the  allcwed 
amount  of  time  and  the  wind  is  not  providing  enough  energy  for 
the  system,  the  farn  tractor  is  attached  to  the  system  as  a 
back  up.  Running  t)\e  tractor  cost  additional  money.  Therefore, 
the  total  amount  of  tractor  enei’gy  used  in  the  complete  simulation 
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1 "U  l:'>  rocv'.v'.oil  in  Livj  lolr.l  l.i-.ictcr  u:;;.!  TTl-U  accuivinlator. 

Ail  " . Ai'-iinY ".3"  ai'C  i;r-,'uao<l  i ii  '< i lu.va  t L-houra  (Ic'.Vli)  and  all 
■’ 'S"  .>  cicj'i  t';;.;.- d ir.  k iivAv.i  t la  (X'V), 


raoluri:  ■."•■i  • o.-r; , ’rhc: 


1 tiona  l!i,iC  dear.i'ib'i  Oflioa 


I a.nd  11.-'  aaai.T.ptiona  nr.derlyif^^  ti.oin  arc  aKplalnad  bolow. 


n'.lr.'.in  tho  fiy.vheel  power  ir  .J  air.cti  tlio  laal  ('■'•ied  plus  tha 
1 (MO  lor  pov/ar  rood  rirK  ? ’ lio  last  pt.'''uc.i,  .■■'.ciinr.i  oai  onoi'py 
and  ovt  oP  tl'.e  -.'.rac  ca:  bi.'-cr,./:'  a roar  or 

Iranatrission  wl  i h an  a::soo '.atr’d  o f f io  iar.v'y . raoro  or.oriiii.a  aro 
'lino  affcolod  'o./  tho  oiiMo irnoy  .-.u'  iho  iy. I \t-o'ol:’.  Ir.'tl.'.l 
valu.'s  nuol  bo  i'rovidLd  I'cr  ail  lovols.  In  tho  case  cC  t'A.i  it  la 
t • oi.i-.ja  that  ti'.o  initial  It-vol  ic  ''i;S  ol'  tha  n.y..inol’r.  onor.sy 
capacity.  A factor  o.'  is  inclucL  i to  convert  tho  ti;r.c  units 
fro;,!  days  to  lic’i-ra. 
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¥ ) 


If  ix  v.,;ro  p'insLble  to  rttcro  an  irfirilto  .■inoaot  of  energy 
in  the  fly  vlieol,  fjy/.'lu  el  L-nci’gy  KiVK  ati>i  the  lot.il  a'/.ii  l.'ibio 

energy  T.VE  ’.voulil  alAayr.  bo  cijofil.  But  tl'.e  flywheel  C-ia  Oiily 
more  n finite  nnount  of  energy.  This  j.TCunt  of  i-tor.ago  in 
i-eferro^l  to  as  the  flywiieeV  energy  e.'iMclty  FV.'KC  and  is  c'Xgressod 
in  Itii'li's.  So,  in  order  to  nssumc  a r.iiKi  ir.u.T,  oner-y  storage  caraoity, 
t’WK  is  tl'iO  sair.e  n:>  '\'kE  n.'.til  i'Aa  b.ioop'.es  larger  linn  iaVKf.  Wh.en- 
over  TA’t  is  ('ro.ator  tirir.  , the  flywheel  is  ‘sr.uwiud  to  bo 

filled  with  energy  to  its  o.'ifaci  ty . 

A f ',)S- .K  MI  nte  WI  0 , r Af-'O  •> 

c fWtC-lOOO 

fu'  - ftY.n(rri.  f'lt.njY  iK.ri) 

giv)  - O'! ru.'i'riiN  I't'-!  ’t'lNivu.’" 

..  F'  Y'-'HL-til.  .NH'-v'-Y  d.'.i’VCiry  IK, .'Ml 

i.M  --  iijr-\i.  Av. maalh  f'li;  v:.y  K'./ni 


3 i'CiVEil  TO  THE  y..r'.vh.'-’hl,  l-Tr'.V 

The  power  to  the  flyv.-liool  i'Tl'W  is  .a  r.ate  that  controls  th.o 
wind  ener/ty  inpe. t to  the  flywiieel.  For  the  next  period  in  tune 
Pl'I-tV  is  assumed  to  be  equal  to  the  r.w'r  frem  th.o  rotor  ri'il 
durin.'i  the  present  reriod  if  IFH  in  positive.  If  rtii  is  neg.ativo, 
Fl’l'W  is  sot  to  scro. 

(}  iMFw.XL  Cl  I'Ml''  '■  .K,  I 

i>i  f,,;  ..  i n.s't  t ui  FLY.''H‘;ri.  i'<i) 

- nYNA'U'  (I’.i'.ir.  nnt.'ellM'l 

(Sit.  I no  I < 1) 

pp:;  _ en.jft  r.'i'M  miTot  iksi 


4 FOWKR  KKOM  Thl'l  RO'I'OR  TFR 

Wind  power,  1’^,  i.s  calculated  frot,)  the  kinetic  energy,  KK, 
per  unit  time,  t. 


whf.o  m i:;  mnas 


.liiil  1*1  il‘i  t'lofj  rate,  V in  viind  spj..'.!,  ; is  nir 
der.-iity,  .iri'.i  A is  cro;jr>-s<;cti  o.i.il  area  parp  ir/l  icvi  iar  to  tho 
flow.  For  the  exanuiiitlr/n  of  turbiiv;  per  foi  niar.co , Lher.o 

fie  tore  eoinbii.e  to  ,‘;ivo  a function  of  the  fiH'r.ii 


= '^.3!j5  ■<  lO"*’  ;;i .Ic  va  ti  s/squ.ii  e ineUer 

wb- I'O  V is  in  Fnol.r,.  Hn./over,  this  ns-riinos  tiiat  the  x’otor  is 
lOO;.'  efficient  and  that  it  wii.l  perform  over  an  infinite  ran,:je 


of  v.-inxi  speeds.  A more  rowilistio  fori  hei 

•'S  ?. 

P (kW)  -•  f.3.S5  X to  ^ :<.  rotor  sise(i7.  ) x rotor 

effioieiKiy  (deciTiai  fraction)  x (V"^  - U ') 


'••.hero  V^,  (in  kr.o  s)  i^3  to.e  n’ind  sp.ic.i  at  '..‘hieh  t,h-:  rotor  wifi 


ovoocome  friction  and  'oe.-jin  to  turn,  is  sore  times  referrovi  ta 


as  the  cut-in  v.’lnd  spotu’..  ! .1  aiditloa  thoro  is  a rated  v.ind 


st'oed  V.j  of  the  system,  beyond  which  no  increase  in  power  will 

r\ 

result  due  to  an  inorease  in  wind  speed.  T'ae  fell  theorem,  (ff) 
indic.rtes  tlia  t Lho  rotor  efficiency  lias  a theoretical  r.axinuun 
of  l6/'.?7ths  (=0.5>'6),  yeasui'od  values  ranye  from  0.‘'t5  to 

c.to  (15). 

Af'.otlier  facto;'  that  s)'.ou!d  be  conaidered  when  evam.ir.itv'  tlio 
power  An  tlie  wind  is  tiie  licif,b.t  of  tlic  windaiili,  Near  the  grour.d 
the  vxind  is  sloi/ed  by  friction.  The  amount  of  this  fiictlon 
dopeiids  on  the  rouithness  of  the  surface  and  obstacle.;  in  t.he  wind 
path.  Studies  have  shown  that  a powtir  law  increas;  of  speed 
with  hoir.ht  is  adequately  d.  joriptivo  for  practical  purposes  lid) 
Over  flat  terrain  similar  to  Ih.at  which  we  are  con- idcrir.t  for 
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iri'i 'ja,  thi'-:  ^o'.vjc  rjLatLof  oan  b..*  u <{  rt-’yaud  am 

f>  z V? 

L.  =z  {?-) 

o o 

VihclO  P iij  the  po'.-er  avail  able,  z is  iijt;  a;i?;lic  above  ground, 
and  ths  subscripted  value.',  ar.}  thoua  at  the  anc'uQn.e  ce r heighc 
at  which  the  wind  (lata  'waa  gr.therc-d.  Thur. , ac  the  I'-oigtit  of  tha 
wif.d  turbine  is  increacod,  tno  power  avillablc  I'rO't  the  wind 
incrao.aoa. 

All  of  thscc  principles  are  cc.cridc-i od  in  the  oquatiens  for 
power  .fron  the  rotor. 
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5 ADJUSTED  IVIllO  STEED  ’ WSA 

As  has  already  boon  mentioned,  there  exists  a wir.-d  speed  above 
wlileli  an  increase  in  wind  speed  will  not  result  in  an  inorcarc 
in  power  from  tlic  rotor.  This  is  referred  to  as  the  rated  wind 
speed  K'.v3  of  the  system.  Tlic  adjusted  wind  speed  '.‘/SA  is  actual 
wind  data  that  has  been  adjusted  so  that  any  wii:d  data  that  is 
recorded  as  being  liighcr  tti  in  R'.VS  is  set  equal  to  RWS . This 
adjustment  is  made  here  in  a separate  auxiliai’y  equation  in 
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r;j-J<.T  to  .'jiiijunf/  the  eqoatioa  for  po'.voi’  '.ro:;i  th.o  rotor 


\ '..r,;,  .■:<  = L 1.  UM  t-'-S  , rtS  . t S . K.  5 

: Rws  ‘>.1 

- >iNi)S'’nii  (M'.ii;sr:n  (■<\";rs> 
Cilr  “ i<  i,ic  1 •J  JCi  ( 0.) 

Ko’S  - WIM) 

rtb  “ bf'ioO  r . I t vsoUS) 


(}  Wr.ND  SPrl'iU  VS 

The  v.iiiT  . ,:eoi  ’a:,-. 
Vy  the  ’.V  th . r 

Vojthor  huro ’•’>  '.vlivi  o >. t:i 
.■'.vrro;:je  ?.s  c'e.; .'•v.v.’.  or.  o. 


■ for  ih.lt;  •prt.fjrai;'.  \:i\3  rc.'orde’l 
la'oau  i:j  Citv,  ihtr.aar.  lorv;  tcxan 

ir;  reoor.;od  as  .r.  o.io  i'!lr..ite  vlou.i.l 
v.-lnd  sree:l  ir.etor.  Tdiis  is  tyj>i •Ml.Iy 


V ao.-hri  every  t.’n'oo  locur.-;. 
1.5  10  .;..vtevs.  ut..ov*:r,  li.o 
irct.rs  high. 


The  rt.^uoari  height  for  r’oerdirq^ 

1 L'->d-jc-  City  Is  only  6.6 


The  ci'.blc  rc.sfoi.so  !<>  oani  :;i'ior-ii  of  the  py.vor  i.'i  the  v.inci 


.’n.a’cos  prsdic  tion.s  very  ocj'.ai  tlve  to  av'.'re.yifvj. 
or-.or~y  cutp’.it  of  a systo  i '.vLll  alv.-'ys  I'o  high.cr 
predicted  ly  nsing  average  wlrAi  spceiis  (l6),  Ti 
can  bo  reduced  by  naln.;;  as  short  an  avcr.agin,;  I 
possible.  The  'rfcathor  hnro.'ui's  3 h.oar  interv.il 
The  wir.d  data  is  represented  i.i  table  form 


The  actual 
(.ban  the  output 
his  iir-crc;,'o:'.cy 
ime  period  as 
is  used  here. 


•vr. rit.M,  -..'iNi)  i'.\T  ‘ j'jn,  jut , anc  trjh 
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I'V/ru 


The  MOo)iainca  I irt'i.j'n  tion  pu:<.p  di'av;."  its  needed  power  from 
tfie  riywJir.el.  The  i'lyvy.hof-l  power  use, 30  FWPU  will  he  equal  to 
lUo  puT.p  power  der.aTvd  PPO  as  lore;  as  there  is  o.mrqy  availablo 
in  the  flywheel.  V.'iien  tlu>  enot-yy  in  the  fJ yvvhce  1 . drops  lelcv/  a 
sr''ciried  low  e;-.er^y  level  I.ET..,  tiio  pr.r.p  will  stop  drnwiitf  power. 
Viher.evi.'r  the  pui.ip  is  siurt  off,  ‘'■iowr."  time  is  aceiuiiulatirij. 

The  .•q.'ste.'n  can  onl^'  tolerate  1 certain  ae:ouni,  of  down  time, 
dc'.iti  ti-.io  allowed  per  A-i]/  OTA.  If  this  amount  cf  dovni  time  has 
already  accumulated  and  the  flj'wheel  energy  is  still  low,  the 
tranter  tack  up  will  be  sv.-itohed  on  an.d  the  puuip  v/ill  begin  to 
draw  rower  from  the  flywheel  again. 

Tiius,  lt»  do  r07.  .^lining  F'fPU  for  ti'io  next  interval  of  time,  cne 
first  oxaminer  the  seutten  to  de  * '.n-.-iine  i.f  the  sys  cem  was  acc".r.vala  ting 
down  time  in  the  protious  ti;i.i!  interval  by  looking  at  the  down 
time  .accumulat  io;'  rule  DTAM . If  no  down  time  was  accumulating , 
then  [■’'APU  will  equal  PPD.  If  down  time  was  accumulating  and  UTA 
has  not  been  exoueded,  tiie  systora  will  shut  down  for  the  next 
time  j.ntorval.  Ar.d  if  down  tin.o  was  accumulating  and  DTA  has 
been  exceeded,  FW'p'J  will  equal  PPD. 
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'I  TKAC.M’OK  1\)\V.:A  IK'.ACK  TI'LI 

'1  I'i'i' I I'l'  iv^W'M'  .v  t'lOK  >U'  fur  ! 1 i * ''U  -v'l 

v.OiiMi  111'’  iviiiii  1:;  Mowin';,  ",'hin  t'.u'K  nti  in  no  t l,vo  ti.wi  when 

l.lu'  I'lywhoi'l  i’iv.'n;;y  I -vnl  I'V.'l':  In  tow  .ni>l  r.yiiU'm  h.;;:  ilrtMily 
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.'iny  I'Oint.  tn  t.  i m ■ , tr.K'tof  ['..'wi'f  in.iy.e  i:*  I'ljo.i'  lo  .•.no  or  tlio 
t ivu;  l.»'f  {'iwiM'  intv'  tin'  t"l  y .'.li'' n I 'I't't  l''’,V , 
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to  ’I'K/vd'o.^  I into  th;':  y'iYwiihj-:^  Ti-'ci'.v 

I'i'iic  loi-  pi.'Wi}r  iato  \;)io  x'lywliool  ia  clcfinod  as  tha  [.'toOi'ct 
of  Itu?  txMotor  ; i.vcf  I'utir.;  'I’  . i!  ami  the  farac  tov-to-f.ly  vhou  I 
ti'anrjmisiiion  off laioncy.  I'at'.i  tractot-s  have  booo.-.'.a  more  powor-ful 
a-Acti  your  as  tho  ot; uipi.ient  tl’.xt  ttui  fari.'.ei’  uses  has  .^ro'vii  in  si/.c 
and  pcv;er  reieiremoixta.  'ro'',ay,  tr.'ic  t;ors  'c/ith  tOO  Lo  200  HP 
ra't  lr  ;':  are  ret  nr.noir.i.ion.  A trac  :a!r  ratin.^  of  100  '/.W  (appro;;!  ...i  toly 
HP)  if,  aesuJi’.ed  for  liiLa.  oyol;Oin. 


ft 

C 

C 


( 1 t-w,  ^ =T I r 


to 


r!>A  -•  too 

t I 

TPIFA 
■I'o;! 
TTf  Hi- 


I 1.  1 
1 0 . 


r^ftCioe,  i-'OaLx  Tir.)  fLY.'n'fL  (k,-;) 

I \ ' r lu  ^0  *1 ; ■<  " A J i 'is  ( X w ) 

rXAO  rOK-  ! f-fL  vrr  M.  UAttMASl  '■'l  e.m- ICtb-lCV 
( !) ! 0 ! Al.  I'r'ACriOM) 


11  TOTAL  TRAC.''0H  '-H.’LiiGY  H'.H'D  O.'THU 

Total  tractor  oniJrpy  ur;od  T'-ili:'.)  is  a level.  Kvery  hour  that 
the  tractor  is  opei  :i  tirri  a’ld  ;'a  t titi.'-;  energy  into  l}.o  irrigation 
syslon  costs  Rorcy.  It  is,  tacrofore,  necessary  to  keep  a runnirg 
total  of  hovf  niuc!\  time  the  tractor  is  curned  on  in  order  to 
cxanlno  the  econoRic,  aspects  of  the  irrigation  systcrii.  At  any 
point  in  time  the  total  trac  Lor  energy  vised  is  e-Tvial  t.o  the 
total  tractor  energy  used  up  to  the  previous  period  in  time  plus 
the  tractor  energy  used  in  tho  previous  period  of  time.  One 
must  pay  for  all  tho  energy  put  out  by  tlvo  tractor,  I'.ot  just  . 
that  energy  that  is  put  into  tho  flyv.liooi  after  passing  through 
tho  trac tor-to-f ly.vlieel  transnission.  Th.erofore,  t'ne  tractcr-to- 
flywhecl  transmission  oi'f ic: i ency  must  be  divided  out  of  t)ie 
tractor  power  used  in  live  previous  tim.c  interval,  'devo  initial 
conditions  are  assumcvl. 
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I.  1 Ttu-K'-TTfu.  .i » t i:r ) ( ri'i.i.,  j(\ ) ".’'t  / r irwii  ii 

N TlrU  0 l.'itllM-  C.l'Ji'I  r M.NS  M.l 

TTl.U  ~ lUT-M.  ittACd)^  gN‘.id.»'  US, 'ID  (.sJ!l) 

DT  " TiM,': 

I ru  - u '.r.r 

TVi'wC:  - r KAi:i  n l-r  c - I L V'Wt! , i 1.  I S.AMOn  ! lOU 

l;'mcicncy  ( cr.c, if'.Ai.  rnAcsniMi 


DOWN  TIV.S  DOWN 

Dov.a  liwc  DC'W.l  Is  a l.ovo'!,  arid  lu  tho  aTnouat  of  time,  i.n 
hour?.,  tiiat  the  i rri t'.  o.'i  r.y:?teii!  v/imlO.  tr  do'.,’n  day,  one 

day  at  a time,  if  there  were  r,o  D.aek  up.  DOV/ji  at  : wip'  point  i.n 
t.ir.o  iu  equal  to  tl'.e  rimount  ei'  down  tin.e  .locu.'i'.ula tod  cince-  tfic 
pvevio.js  a idni^jliv  piii.i  liie  .aieo^.nt  of  dowr'.  time  -'.li.'.t  wi  ti  ooour 
ii\  ti'.o  neat  in',  etv.il  . If  tl; ; fJy.vliecl  e.iei-^^y  level  lov.-l  iu 
holcw  t)ie  cpeoifi'.'i  low  cr.ery.v  level  I.lil,  at  a point,  in  tiiae, 
tf!  .-n  : i.!itic:r,il  down  time  will  be  aec 'M',- i.-', tc.i  duriy,'  the  ne-tt 
tviae  , interval.  .\t  t!io  -'nd  of  ti;o  day,  the  dov.ui  time  o.ccur.nlatiu' 


is  r;et  equal  to  zero,  'I'hia 
entire  function  for  DO'WN  by 
do'.'/n  time  is  assi;  r.(!d  .>.,5  the 


if  .ac'-'o.-.pl  i.'lied  by  multi  p^yj.i’^^  tlie 
a '.yeio.a  nv.atu.i  index  SDI.  Zero 
initial  v.il'io  of  DOWN. 
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13  DtWN  Tn.;/;  Af.Oi.ir.ui.A'i'iCN  ka'I'd  otar 

Tho  down  time  .acouMul.a  lion  rate  Dl’AR  is  an  indc.x  of  v.dietnor 
tho  jiystem  is  accunul.iti  (J.'-'wn  time  or  not.  If  tiio  flywheel 
enei'f.y  level  FWK  is  below  tho  specified  low  onerf,y  level  LFL, 


lUon  tlov.-.i  wMl  ho  aoouri'u.l.i toJ  nt  a x'ato  rolatcl  to  tho  pomj) 

pov/oi'  PFD.  liTA.'  ilooa  not  dvroctly  inJ'lucnoo  t)'.o  lovcl 

of  Down.  fi'.ntead,  it  as  an  index  aa  to  tho  immediate 

future  accumulation  of  dov/n  time. 
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I'i  SYhT'd.Vi  STATUo  Ih'Di-JX  h'Sr 

The  aystcir.  atatua  index  [ resota  ti'.o  00', '/o'  aocu.'iui.Iator  to 
7.oro  every  24  ,‘iourt;  .t  iiidni.jht.  It  in  alno  uned  in  tho  calculation 
of  the  total  do'vn  tn,-..-;  accM...\iIation  rate  'CDTAR . It  1.30  the  value 
of  1 at  every  point  in  time  except  at  i.iivinigti !.  each  day  sviien 
it  )ian  tho  value  0.  Tixenc  value;;  nve  n Lured  in  table  form. 

Figure  3,2,4  shown  on'-'  section  of  this  t-.Me  -.s  the  computer 
interprets  the  system. 

A SSt.K^TAntri  SSI  r . U‘,)  l-'t 

T SS1T;0/ 

X i / I / I / I / W1  / W J.'l  / I / I / 1 / I / I / I / J /I  / t / ! / 1 / I / l/I  / J/ 

X 1/  1/  1/  1/  5 / 1/  1/ J/1.  / I / l / l / l / I / I -'J/l  / I /■  I / l/l  / I /I  / )/ 


X 1/  1/  1/  1/  ! / 1/  ! / O/  1 / ’.  / 1/  1/  i ' 1/1  / J / 1/  1/  1/  1/  1/  I''  !/->/ 

X 1 / I / 1 / 1 / I / ’,  / 1 / 0 / I / I / I 1 / 1 / 1 / 1 / .)/  1 / I / 1 / I / I ' V / 1/  d 

SSI  - svsnv  STMUS  M’;\  ( t;l 'NiS  ir\L  fSS  ) 

SSIT  - SYSli''  SlATtS  I'ileX  r.ihtl 

t5  ENF.RrjY  SURPLbo  ('.‘.•IND)  Ir.t  INDLX  KSWl 

In  a given  interv.al  in  time  it  is  very  possible  that  tho 
enei'gy  available  from  tho  wind  exceeds  the  available  storage 
remainiryi  in  the  flywheel.  In  such  a cane  the  excess  energy 


can  not  bi*  pvit  to  vr^o. 


Thi:j 


uri'lus  v/l!vl  onoi'cy  ia  noted  by  tbo 


anxiliary  ro.^iat.T  hbWl.  It  in  nimply  any  poaitivo  d U'roi-onco 
between  the  tct.il  av;niabln  oroc'i'^y  TAH  and  th.o  flywhool  .niorf^y 
oapacity  FWi:0.  Thun,  ones  tins  f* yv/hoel  is  liiled  to  capacity, 

■".'iy  excess  wind  enoiy'y  in.  noted. 

\ : Swi  I -tvi-r.  ,0  I 13 

(-s.ii  !S  U5LP  ro  I VAi.  b'.  C'-  fHii  I'.nrot  31/.F  - 

riY^tn-a.  cvnuNviioM 

I ', ,11  - .-tn.’id/t  .tu.'i’i'.’i  .M'jn  (isi  ri'.'f  0 (.twii) 

y/i.<  - ov'i.-.po  l UMC  fifp  r-dt  >M.\atdM' 

TAc:  - I or  At  AV  !l. /'.''I.  c '"n-,  t’.Y  (KOO 

rw'tr.  - i i.Ywiirti.  t.’jaa.Y  CAntciiY  ia.mM) 

16  TOTAL  J.'.Id'/adY  SURPnlT.  (b'Tiai)  1st  .INnitX  TKSWl 

Total  coi?7'gy  sunrtlur,  TL'.T.s'l  acc'.uinilates  a runniivT  total  of 
FSV.'I  while  the  systor.  in  oponcti nj.  T)i5.s  lovol  is  equal  to  Tliv/.il 
up  to  the  previous  poiiU  in  time  plus  the  surplus  wind  cneiqqy 
(tst  index)  of  the  previous  period.  The  units  are  kiicwat t-hours 
and  ccro  initial  oonditioriS  are  assumed. 

By  examining;  the  total  energy  surplus  (wind)  1st  Ir.Jex  and 
the  total  tractor  energy  used  TTK'U  .\t  the  cn.d  of  a simulated 
stiir.ner,  a logical  decision  can  be  made  concerning  possible  raod- 
ifieations  in  the  systoin.  If  TKLV.  1 and  TTbtl  ax'e  about  equil,  one 
would  know  t!\ac  over  tlie  summer  there  is  onoudii  powt.-r  in  the  wind 
to  moot  the  pov/er  doma.ads  of  the  ii'rigation  system  but  that 
possibly  a larger  flywlieel  is  needed  to  even  out  periods  of 
low  and  periods  of  high  wind.  Likewise,  if  TTEU  is  nuch  larger 
than  TKbWl  , tficn  a larger  rotor  is  needed  to  make  tlie  syst-nii  loss 
depor.dont  on  back  \ip  energy. 
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I f "iwi  .N --1  i s-v  1 .,1  ^ f ''r ) 1 , I ) ■'  i 

ri'jwi^o  riiri  ii  co'nniiMS  i 

Ti-S'-ll  ' liiTAi.  i:,'j(Pr,Y  SU'-'LllS  .-1'''.)  tlST  ( ••  J t<  I 

or  - ! iM(  nrsovV.M.  ioayS) 

:sKi  - sut:''i.us  .frio  iisr  noi;x)  kvoi 


1?  KNl'RUY  SrUMM.'Jj  ^nd  ICT'KX  E'3’^2 

\ wlr.d  pow„i'  ':;y::toni  oi'  aj.y  kiod  •/ill  .r  'Alr.d  spoad  at 

•.vhioli  thir  ayatiin  '.vili  oporato  at  i ta  p-vak.  Thii;  rpoed  la 

rorao'cad  to  aa  Cba  rated  r/ind  apaod  aiirca  thia  ia  U-.o  'j-pood 
at  v/hich  the  t-y.-j  p,  rforai  at  ita  .■.axiinam  ratir?’:.  'C'to 

energy  available:  Irora  tl'o  wiad  ia  diroc'cly  propvrtlor-ril  to  the 

cube  of  tiia  v/ii'.l  apood.  Thna , aa  the  v/i;vi  opeec  ii'.crea'aeo , i-.ove 

a;.d  n-.oi'c  eaar,;y  can  be  or.ti'acted  by  tive  -vird  pc-.er  ay.3to.-n.  Hut, 
na  has  boon  rxpl.xlred  alro-uiy,  beyci.i  the  rated  •■.■.'nd  opc-'d  of  lire 
nyaton,  no  n-l-.l i I \ or.  ' I onei'py  ea'a  bo  derived  fror.  an  j.r.:r.'';-..-30  in 
w.ind  speed.  This  2r.d  iiidex  of  surpl'us  ’.vlnd  energy  koepr.  track 
of  the  amount  of  oi. ir.  fao  v.'ind  t'l.a  t Cuur  .not  re  e.xtractod 

due  to  t)io  rated  '.vir.d  ;5p‘'ecd  of  the  aystem. 


r S If  .< 

■=C,l  1 ’ 

t C ; 

1.  > •>- 

V I-;  E 

( ' .1  • K.  • ' 3 • X r n > . K ' ■ ( 

. . ? :W  ■ L 

lb  1 

-tlSs., 

1 C 1 

> 1 ~r>r 

^ , 

. 

•■1  S . K 

, Kvb)  11 

c;  S •< 

1 s 

1 ii-ti 

T'' 

v'  •*« 

fj.-I 

f I HE  ■v\rE  r Vi 

b E E 

A rr  ni, 

'» V S 

r ! t 

E b V,  r 

- 

i 

SP-'f'Lub  ..I'V)  1 1 

10.  \) 

( <.au 

nir 

' 

I'.y..  •. 

f.-  J 

ten 

n 

- 

1 Tin 

’> 

‘ M' If,  i'roY  (■.)rcr“'.L  ! 

N 1 t 

If  VI 

() . 

3(,r,-S 

- 

C -'lb 

i 

\T  C 

ill  n I ,:r  SA  r m on  , 

vr-  ’ r 

ro  [>t. h- 

.> 

'j\  i i 

KS 

- 

P.OI  1 

p 

S W f 

'Ulfhbl 

'..'S 

- 

V ! N n 

;>!’  1 1. 0 

l'.'.  !.•.  (\ElIt) 

1 xa( . A 

I!  '3  ! 1 

l.e 

, ( HSr 

1 

- 

Hi  I vi-Hi  fM  r.-is  r./vf 

rjo'i 

HSf 

HE  i 

1 r 

tr.Ml'O'.  (.\r.  t K .VI'. 

L"  I l 

l M"  I tllT 

( 

l,'! 

1 

' ^'S  ) 01  Y 1 ')i  DAY  \ 

•;i  V,  ' 

T:A<  1!  I 

■;  i 1 

C I.S 

- 

CU'  • 

1 %’ 

-S-'E 

fi:  (v'jnrsi 

I’t  r, 

- 

I’d  Vl 

n 

f-Ou'I 

iiH  V ir.M  1 ■<■.;) 

i>  WS 

■;  A 1 i-. 

.) 

W1  NO 

.ShtED  1 K J IT  S ) 

18  TOTAL  EWIHin'  tiL’KriUti  (k'rND)  ;:nd  I.NhoX  TfSW2 
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Tc(.<].  energy  .-J'.u  plu:;  TK.";V/2  ocr  i;-.,ti  1;>  Lo  i a ruf'.ui  p;;  tcl  il  of 
LoWS  tho  aimulitiou  of  the  ay:; ten.  Thin  iovol  is 

CTunl  to  T;ioW2  ut)  to  tJ'.'?  prcrioi:.';  point  in  time  plus  Dio  nurpius 
v.-iiul  ci\ocQ'j  (2n;l  index)  of  tho  [u'cvious  poi'iod,  Th)  units  nro 
kilowatt-hours  rnul  /.ero  initial  o.onJltions  are  cissoMoii.  As  )ias 
niready  boon  cxpli.inod , chore  is  a limit  as  to  hov/  lar.'O  tho 
rntod  vind  r.pood  of  a systor,  i.ny  prao  tio."  lly  bo  nena Idored , 
Therefore,  it  may  nove-.'  bo  possible  to  cbnr.^u  HV/b  or.oiiy.h  to 
reduce  TnS'jVS  to  xero. 

( r f >.K  D(  S'.'d ..)  oi ) ( r s c,’ . j ) -^  d it) 

iNi  n:sw2--o  Di.i 

- inrAt  :ytra.Y  SU'^ei  n;  '..I'P)  (S'.S  I'lD^X)  (CV'i) 

i)r  - 1 ru  inrtrv.H  i :)\ys) 

ESVP  • CN.-OoY  S'J’Ji’eJS  JIN.)  (dbO  I'hl'.o  (it'.JD 

19  (TOTAL)  CAPITAL  COST  CG 

Capitol  cost  CC  is  expressed  in  doliar".  CC  'r.sluuca  tbo; 
cost  of  the  flywheel  and  rotor  system.  Iti  aa;iition,  each  time 
tb.e  tractor  is  tirenod  on  as  a back  up,  tho  cost  of  this  enoryy 
is  assessed. 

The  rotor  cost  is  assumed  to  bo  $100  per  square  mstcr  and 
includes  the  cost  of  a 30  meter  lower.  This  fip,urc  is  based  on 
cost  information  provided  by  various  maivifac  tuicrs . It 
assumes  Die  i-otcr  will  bo  a Oarrious  desiryi  wiDi  a Savonious 
s o 1 f - s t a c t i n.t  c a I'.a  b i I i t y , 

Tho  flywheel  eapi lal  cost  is  assumed  to  be  $35  per  kWh  of 
stora.'ic  capacity,  fruch  a flywheel  would  be  approxin.a to ly  G - 10 
feet  ill  diamoter  and  wci.qh  boiwcrn  .50  and  75  tons  for  a stora,T;o 
capacity  up  to  5000  kWh's  (19). 
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The  other  cor:t  conoidercd  that  cf  the  back  up, 
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.'tor  ‘CnerfT,/. 
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of-  . 

.uierfrv  coats 

for 

tiro 

purposes  of 

A:;nu.T.in,3  a craetor  rating  oC  lOO  each 
between  .;'.0,10  and  :i;0,20,  t'0.20  io  anaun.ed 


thlc  v.odel. 


It  Is  not  tlie  uithor'a  intention  to  any  ci'  the  oteor 

i.oe,  ta  aonociatod  v.'i^h  a synten  such  an  in  proponed  here.  It  i.a 
ancnm'.a,  though,  that  t!.e  oont  of  the  trancri.lsnj.onn  and  tlie  ^ear 
box  are  tr. inor  atvi  tiiey  av-a  trioeefore  riot  oon.^idcri.'d . Th.e  rotor 
car i tel  cost  and  r.ha  flyvlieel  capital  cost  can  be  asnancd  to 
equal  7cro  if  Just  too  cent  of  niai ntalri;\'  the  back  up  systein 
is  d‘-‘oLrod, 
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of  . K-T  T t-r . tf  :r,  rf<S”'rCC  rf  'Ibc  I ? 

r f r c = . 2 0 1 . 1 

3cr  = rJd  i‘>.2 
f /<  d C - 3 j t 1 • 3 

cc  - irnrAt)  ffo’iTAi  c.isr  (*,) 
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Tree  - fctcrnK  er.'eKdv  Cif-’iT.M.  r-isr  i t/k'.: 
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FKCC  - f-LYel'ifL  capital  lOS  I (t/itAH) 


20  TOTAL  DOWM  TlliE  ACCUliU I.ATIOti  RATE  TLTAR 

This  rate  controls  the  accuir.ulation  of  the  total  n\:  .;ber  of 
hours  the  system  is  actually  "down"  dur.ii\>3  a simulation  run.  As 
the  systcun  has  been  defined,  there  is  a ir.riximu'’.'  nuirher  of  hours 
of  down  time  r<ll0'.;.''d  per  day  DTA.  At  tin;  e:id  of  each  day,  TOCAR 
looks  at  the  DC'./N  aocumvilator . If  DOiVN  is  less  than  DTA,  TDTAR 


assunies  the  value  ef  DOIVK.  If  DOV/h  is  , greater  than  DTA,  TDTAR 
assumes  the  value  of  DTA.  In  the  calculation  of  the  total  down 


tiT.c  DN'TOT,  TDTAR  is  multiplied  by  the  time  interval  DT.  DT  is 
defined  as  borne  equal  to  .125  days,  or  l/6th  of  a day. 
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2t  TOTAL  DO'.ViN  T uAE  ACCIi^.iULATE'O  DNTOT 

This  lovcl  accuip.ula tc-G  cho  total  n'j.T.bar  of  hours  the  system 
is  "down"  during;  a simulntion  run,  Di-iTOT  can  be  used  to  determine 
the  load  factor  (percent  cf  time  the  system  is  oneratiit;; ) dui'ing 
a Gimuiation  run.  Zcr.'o  initial  co,.)’. tionn  ar-i  assumed. 


L Of.  ••'’T  .r. -o:.|T '!T  , J t ( d!"  I M O ! A'^.  -i'. » 21 

S OiTorrj  iMTl.n.  Ci'hor  I r.s  2l,l 

ON'TO,'-  - total  OC/iN  '.CC  J ojl.aT  >20  (rP'ij'-'SI 

01  - rite  L.MtRVAL  I'TA'Tt) 
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22  LOAD  FACTOR  LOADFAG 

The  load  factor,  as  oalculatc.l  by  this  auxiliary  equation, 
is  the  percent  of  time  the  irrigation  system  i.s  cpcratir.,T  durin/; 
the  summer  months  of  Juno,  July,  and  August.  At  arAr  point  in 
tine,  LCAOFAC  is  the  percent  of  the  elapsed  time  that  the  system 
has  been  operating. 

Usin,T  a pump  rated  at  100  kW,  it  would  be  necessary  to 
operate  the  puj/ip  50  “ of  the  time  in  order  to  provide  6 to  8 

acrc-inches  of  irrigation  per  month.  The  down  time  allo'.ved 
per  day  DTA  can  be  adjusted  until,  at  the  end  of  the  cutniner, 
the  ovcr.ill  lord  factor  is  between  ^0  and  60;i, 
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T tit.ro.i'if  tr  icn.  A scliom,^  tie  oT  Option  2 npptjnra  in  3«3.1« 

111  this  f.;ystoM  .ho  wind  tr.rns  a wind  turtino  (vvirvJiiii  1 1 } v.jiich  is 
coupled  to  ai  o'ooti  ical  altornatoe.  The  elec i.i-io i ty  ^oi'.ci-atod 
is  then  reclitiod  and  I'rd.  into  an  el  ;o  trolysi  s coll.  The 
electrolysis  coll  prodi-.eos  ''.yaro{toi\  and  oxygen  Troni  v/ator  (11,0). 
Tlic  h\drc|;cn  gin  predueed  is  then  pie:d  off  ruvi  stored  in  an 
K :chae,;!  ted  natural  gas  well.  V.'lien  it  Iccornos  r.ecene.-'.ry  to  provide 
er.oiy-y  for  Lrvi  ;aticn,  the  i'.ydrogen  is  piped,  to  an  Intern.al 
conbuf.  ticn  on^jine  wluch  drives  the  irrip:ation  punp.  Pronar.e 
provides  a back  up  for  operation  of  ties  system  \.h!;n  r.o  jiydro._,.fn 
is  available. 

Option  2 involves  throe  energy  cor.vurs  ions , wind  powier  to 
electrical,  electrical  to  production  of  hydros^n,  and  burnin)? 
the  hydrogen  to  provide  nicci'.ar.ical  energy,  co::’.p,\red  t.o  the  simple 
conversion  of  wir.d  power  irito  meehaaical  pow’^-r  as  in  Option  1, 
Tims,  Option  2 is  ir.hercntly  loss  efficient,  liov/evcr.  Option  2 
provides  a means  I'or  storing  much  more  energy  at  a lower  cost 
than  Option  1.  Therefore,  the  system  could  easily  operate 
d\iriTvt  tlie  ojitiro  icar,  storiiv:  the  h.'.drogen  pr. -luocd  during 
the  month.s  of  September  through  Wny  for  use  during  June,  July, 
and  August.  Tlui.s,  it  is  possible  to  provide  the  same  amount  of 


total  energy  as  Option  1 at  a lower  total  cost 


electrical  alternator 


IBUSTION 


hematic  of  Option  2 


9) 


A cao'ual  loop  (iia;jr.i.ii  of  Option  2 in.  sl'.o  vn  in  V't^parc  3*  3*2 
and  the  i'i'IwV.-’O  flow  dianjrars  for  tho  model  io  sh  <.vn  in  Figure 

3.3.3. 

§..Z? '^05L  ' Thore  are  four  loveln  i.n  Option  2,  Hydnof^en 

storc-.^c  level  U^SL  is  the  Jimiiber  of  kilc.va tt-houi:.;  of  hydrOji^on 
enoryyy  that  are  in  the  natural  gas  .veli.  at  an.y  given  point  in 
time . 

'j'otal  cnoray  surpius  (vvir.d)  2nd  Lr.c'.cx  'ri'aiv;?.  serves  ti-.e  same 
fnnetior.  as  ThtjV/2  in  Cptior.  1,  Similarly,  dov.n  time  COi/iJ  serves 
the  t-3.nie  fun,tion  as  DOWW  In  Option  1. 

The  prop'ane  used  PU  is  similar  to  the  total  tr.actor  energy 
used  in  Option  t,  V/hen  tlia  system  h-as  boon  "d.cv.ri"  for  the  aliovvcd 
ai.iOi.int  cf  time  and  thorn  is  not  enough  hydrogen  to  power  the 
syst.om,  propane  is  used  to  power  the  irrig.atlon  p'lnp.  F'J  is 
t'le  total  amount  of  propme  usod  up  to  any  point  iii  time. 

Due  to  the  .large  capaoity  fca'  entrgy  storage  cf  Option  2, 
it  io  no  longer  ori tical  to  be  pumping  vatcr  whenever  the  wind  is 
tlowirg.  Therefore,  the  appro.ach  to  cs  cablisnis-.g  the  load  factor 
of  tiie  system  is  different  th;),n  in  Option  1.  Those  differences 
are  explained  below. 

All  "ENTiP.Gf ' d"  .arc  orpro-ssed  in  kilow.att-hourrf  (k-Vh)  and 
ail  "POiv'k'tl 'S"  arc  expressed  in  kilowatts  (kvV). 

Structure  and  .i.-jr.:  oti  onr. . The  iS  equations  that  describe  Option 
2 and  their  underlying  assumptions  are  explained  below, 

1 HYDSOGEN  STOR.IGE  LEVEL  H2SL 

The  hydrogen  storage  level  H2SL  contains  the  number  of 
kilowatt-hours  of  hydrogen  tiiat  are  available  for  use  in  the 
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F'  Aure  '’,^.2  - Casual  Looo  Diap:ram  of  Option  2 
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Figure  3»3.3  - Flow  Diagram  of 
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hytii'0{;eii  intei'naJ.  ccm’ourj  ti.on  on"lno,  lll’SL  a t ai'.y  poLit  i a ti  .".tj 
3.S  r-qual  to  the  value  at  the-  !)reyiou:'j  point  in  time  plus 

the  hydrogen  produced  durin^j  the  last  time  interval  minus  the 
hyiirotpen  used  during  t’ne  last  time  interval.  Tn  addi. tion,  any 
propane  used  as  lack  up  fuel  must  be  considered.  The  amount  of 
hydrogen  into  and  out  of  the  ijas  'well  storage  is  also  effected  by 
the  hydrogen  storage  efficiency  Hk‘JEFF,  The  calculated  value  of 
at  each  point  in  time  is  first  ccinoared  witti  zero  and  the 
larger  of  the  xv/o  is  assigned  as  the  actual  iiSSL,  This  is  to 
ii'.suie  that  there  is  never  any  "r.egativo  energy"  stored  jn  the 
gas  ’.vall.  Zero  initial  constitions  are  assumed  for  this  level. 
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2 HYDSOOMN  PRODUCTICW  RATE  1[2PR 

The  rate  of  production  of  hydrogen  over  an  interval  of  time 
is  assumed  to  be  equal  to  the  amount  of  hydrogen  produced  by  the 
electrolysis  coll  during  the  same  time  period, 

a H2i’^.Ki,=h.?rc.k  2 

H?pr(  - HYOU'IGt^  F.tCOUC  I!  f V TATF  (K'.-!) 

- HYO'-UGKV  FFf'  1 n.  ( tr  ^u'LYb  IS  ) CtLt  I K d ( 


HYDROGEN  IROM  (ELECTROLYSIS ) CELL  H2KC 
The  amount  of  hydrogen  produced  by  the  electrolysis  coll 
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is  equal  to  the  :;.,iovint  oC  eloo  trical  pc\/oc  proiiueoit  'ey  the 
rotor  EPh'il  times  the  rectifier  cri'iciency  Kl-FF  (the  rectifier 
converts  tlie  altern?.  ti  tu;  current  pcv;er  produced  by  the  alternator 
to  direct  current  pov;ur  for  t/ie  production  cf  hydrogen)  time-fi 
the  electrolysis  cell  cfficienoy  C£FF.  The  .rectifier  is  a 
solid  state  device  v/ith  an  efficiency  ea.s.ily  o.uproachiit'  9l5/J. 
Tiierc  presently  docs  not  exist  an  olec  troJ.ysi.s  coll  of  the  size 


needed  for  this  systen.  Ho'vcvor,  f;ener..il  Klectric  projects  a 
cost  of  .jJO/hV/  at  efficiency  for  this  size  of  cell  by  1985 


(20). 

For  this  simulation,  an  efficiency  of  90;1 

is  assumed. 
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The  mechanical  output  from  the  rotor  is  used  to  turn  the 
alternator  to  produce  electrioal  pov/,er.  The  electrical  power 


being  generated  by  the  rotor  iSPFR  will  equal  the  mcclianlcal 
pov/cr  from  the  rotor  times  the  alternator  oi'f.icicncy  AKFF, 


Ef’FR.K-CL  I P t PFR-  K ‘‘l  EP’'  • r '♦ 

EPr-B  - ELsCTRICAI.  PiIvitR  PRCa  ROTOR  (K'Wl 
CLIP  - UYAAMC  LOGIC  RiJNCTtOO 
ppp  _ poser  from  rotor  (K'h) 

AtfF  - ALTCRhAraR  1;!- 1- I C I EN  CY  (OOCIMAl.  rP.ACVlON) 


5 ALTER MATOR  EFFICIENCY  AKFK 

The  alternator  efficiency  is  load  dependent.  A typical 
efficiency-load  curve  for  an  alternator  is  sliown  in  Figure  3.3.*T. 
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In  thin  nyf.Lom  tho  alternit<n'  vvilj. 


ot'oraLin,'  under  full  load 


at  the  rated  wind  npoed  RW^J.  Tiiun , tho  input  ratin.'j  of  tiio 
alteri'.ator  AI.1'RAT  can  be  calculntod  by  computur;  hov/  nuch  [.ower 
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6 roV/RR  RROy,  ThK  RCT.'R  PFR 

FFR  in  Opi.ion  2 is  the  sa:no  as  in  Option  1.  The  effects  of 
rotor  nir.e  and  efficiency,  raced  wind  speed,  tower  height,  and 
cu t-in-rpeed  are  all  ccr.niderod, 
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T'lO  ivind  rp'^e:}  ir::d  in  OpLioi  2 is  the  SDino  ay  ti'.at 

used  in  Option  I v/ith  one  o\co'..tior..  Cptior.  1 is  onlj'  cor.ocnned 
vflth  June,  July,  an.i  A.unu:;',.  But  '.iio  ?.-i3te;r.  ir.^.'iclcd  in  Option 
2 v.ouid  operate  s''-’''r  round  and  .Jmii  uces  wind  ‘.i  Ti.  "t  1 JL*  0 in  all  12 
months.  To  prevent  an  ovevricw  of  the  computer  memory,  wind 
data  is  used  in  3 month  bloeks  ho/sin..infi  with  B'eptemher,  Ooto'cer, 
ttnd  november  and  ending  with  June,  July,  and  Auyust.  Thus,  the 
syisteir.  is  able  to  produce  aiui  store  I'-juirofCcr.  for  9 notilli;:  boCavo 
tlio  irrigation  pr.mp  is  nirnod  on. 
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9 IL/nHOGEM  UGAfJE  RATE  H?-UR 

Ti;c  ):y(!i'C£-!)ri  uca^S”  rato  K2UR  is  sinii;j.x'  to  the  flyviOxueL 
po’.vor  ijsage  Fa'Pu  oi'  Option  1.  K2UR  ’.vi],!  ha  equal  to  the  py.np 
, po-.-.ar  den-.and  FFD  ns  long  as  the  ar.ount  of  hydrogen  in  storage 

' is  greater  than  a specified  lo'.v  civjccy  level  LEI.  Whan  the 

hydrogen  stora.3e  level  HREL  drops  beJow  LKL,  the  puir.p  will  cease 
' bi!;,-ni.’ig  hydrogen.  V/ha.;ever  the  pump  is  turned  off,  "down"  tine 

. is  ncoumulating,  The  ayrrtor,;  can  only  tolerate  a sot  amount  of 

down  time  per  day  DIA.  if  the  cuncunt  of  down  time  accumulated 
is  equal  to  or  great tr  than  DTA , the  system  will  begin  to  burn 
I propane  as  a back  up. 

* Theiefore,  in  do  tor:, lining  }'..?UR  for  tl'.c  nc.<t  intei'val  of 

tine,  one  first  examines  the  systen  to  de  term  I .no  if  the  system 
was  aocumulating  dov.-n  tine  in  the  previous  time  interval.  If 
no  (iown  tire  was  accumulating,  then  !f.?UR  will  equal  FPO,  If 
(lov.n  time  was  accumulating  and  DTA  has  not  been  exceeded,  tfie 
system  will  shut  dov/n  for  the  next  time  interval.  And  if  down 
time  was  accumulating  and  DTA  has  boon  exceeded,  1I2DR  will  cqu.al 
PI'D. 
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10  ru.vp  FOV.liR  DEMAND  PPD 

Tho  pump  power  c’crnarJ  of  Option  2 differs  from  that  of 
Option  1.  In  Option  lit  v.m3  desirous  to  have  thu  system  pumping 
V.'heievor  tho  wind  was  blov/ir^.  Due  to  the  very  large  (assumed 
ir.fiiiite  for  purposes  hero)  storage  capacity  cf  Option  2,  this  is 
no  .longer  important.  Whenever  the  wind  is  blc.ving,  hydrogen  will 
bo  boii'.g  produced.  And  whonover  the  system  is  pumping,  hydrogen 
v;ill  bo  bt-ing  consumed. 

As  has  a.irea.dy  teen  nentionad,  with  a 100  k'Y  pump,  the 
r.yste.m  would  need  to  ho  operating  approximately  bO/^  of  the  time 
(load  factor  = .6o)  in  ordox*  to  moot  ti.o  in  igax  ion.  rou-.ij  roi.ionts. 
Since  it  is  assumed  hero  th.nt  it  makes  r.o  liifferai'.co  when  th.o 

pump  i^5  on,  as  long  as  it  is  on  bO;.!  of  th.o  time,  the  pump  I'cliitg 

i;;  i'cdtioed  to  6o  kW  s.nd  fT.l  is  set  ccju.al  to  r.oro.  So  instead  of 
having  a 100  k'.V  pump  operating  6c;J  cf  thiO  time  as  in  Option  1, 
this  system.  aiJcumos  a 60  kW  pump  operating  100,5  of  the  time, 

A PP3.K  f. L 1 P(  :i  1 ht  Cf  I’fi/T'crr  , r ins.  K , t akoi-f  ) , o , u M'-.k  , l kron)  lo 

PR^60  10. 1 

PtFF=.3  10.2 

lhR(;N'=2f6  10.3 

lRR0FF  = ?'if>  10. A 

PPD  - PDdP  PC'aE.H  CFMA.ND  (KW) 

CLIP  - DYNAMO  L.litIC.  FU.NCTION 

PR  - rOMO  RMIN!'.  (K.YI 

PCFF  - PlIMP  nriCltNCY  tDiClMAL  FKACriUNl 

IRROFF  - Tl^G  IRRIOATnN  IS  TllK.NCD  OFF 

IKRON  - TIMS  IKlUOATItN  IS  T'JPNFO  CN 
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U i'HOt.u  UJAGK  K.Vi’K  PUR 

'I'ho  pr  All'"}  u:;aj;c  x-ato  PUR  coiTcapon la  to  the  tractor  pov;cr 
uaa''o  in  Optii..'ix  1.  Propane  aorves  aa  tlio  back  up  for  Option  2. 
This  tack  up  xa  activated  whoa  the  hydrogen  storage  lovei  ]{2SL 
falls  belov;  a sped  tied  low  oactvy  level  Id-.I,,  Tliercforo , at 
any  point  in  tine,  PUK  is  equal  lo  soro  oi*  PPO. 
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12  PROPA.NIi  USiiO  PU 

Froxiane  used  in  a level  nr, cl  aecut lula tee  ,.t.?  total  anount  of 
propane  used  during  the  s inula tion  run.  b'axxh  kilcvvalt-hour  of 
p’ropar.e  used  cos  fc  money  and  r.ust  thus  be  conddorod  in  the  total 
cost  of  the  systen.  FU  at  a point  in  tine  is  equal  to  the  value 
of  PU  at  the  previous  point  in  tine  plus  tlio  anotnxt  of  propane 
used  during  the  last  Line  interval.  Zero  initial  conditions  are 
tissuined . 
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1 3 DOWN  T [ DOWN 

DOWN  is  defined  the  same  as  in  Option  1,  However,  since 
Option  2 is  designed  to  operate  continuously  during;  June,  July, 
and  August,  no  actual  down  time  will  accumulate.  Tlio  DOWN  level 
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is  kept  in  tlie  ui'Of^rani  to  ir.olio  the  niociel  iiioi'o  flexible  and  to 
have  Option  2 parallel  Option  I, 
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14  uok’N’  Ti;.;tA  lOCunui./iTioi'!  hath  otar 

O'l’AK  io  the  P.ai.AO  in  Option  2 as  in  Option  1.  AGO.in  it 
sorvca  an  an  ir.dex  of  wtetlicr  or  not  the  Gysta.Ti  io  accuKulat j tvi 
down  tiir.e  v/Iiile  rot  actually  directly  : nil  acre  ij'.,"  the  level  of 
DOWN. 
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15  SISTKM  STATUS  INDEX  SSI 

SSI  resets  the  DOk’N  accumulator  to  ccro  every  24  hours  at 
midnif^l'it  the  same  as  in  Option  1.  Eijxnre  3-2.4  shews  one  24  .hour 
pt'riod  as  tiic  computer  interprets  the  table  foi'iiiat. 
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16  EWISRO'X'  3'JRFRLS  (Wi::/))  Rnd  INDKX  ESV/?. 

ESV.’Z  Is  '.he  ~ano  as  in  Option  1,  noting  tlio  a.iiount  of 
crur^Fy  in  tha  winri  that  can  not  bo  oxti’actod  duo  to  tlio  rated 
wind  speed  of  the  system. 

Since  infinite  storat^e  capacity  is  assumed,  there  is  no 
need  for  an  aSWl  as  in  Option  1. 
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1?  TOTAL  lih'ERGY  SURPLUS  (WIKO)  Gnd  INDEX  TKSW2 

THSV.t'!  keeps  a runniipx  total  of  S3Vi'2  through  the  simulation 
I'lj/i  a.e  in  Option  1 . This  level  is  equal  to  TES'.V2  up  to  the 
previous  point  in  time  p.ius  the  surplus  v;ind  cnei'ciy  (2nd  index) 
of  the  previous  period.  The  units  arc  kilowatt-hours  and  zero 
initial  conditions  arc  a.ssuined. 
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18  (TOTAL)  CAPITAL  COr,T  CC 

Capii-;’..l  cost  CC  is  exprccseil  i.n  dollars.  CC  inclu.los  the 
costs  of  the  lotor,  alternator,  electrolysis  cell,  aiKl  the  pa'opane 
used  ajii  propane  tan'c.  The  rotor  capital  cost  is  the  sawo  as 
in  Option  1.  Alterna;  or  coots  in  tl\c  noi,q;hbojraood  of  .f.^O/kW 
rated  output  are  con’.:?on  today  for  alternators  in  the  site  range 
represented  by  ALTRAT.  The  electrolysis  cell  capital  cost 
($bO/k'.V)  is  a projected  figure  from  General  Electric  for  the  year 
1935  (20),  Operating  costs  include  the  cost  of  the  bock  up, 
Prop?u\3  tarJvS  in  the  neigltborhood  of  5OO  gallon  capacity  (15i‘i00 
Id'i’h)  lease  for  around  :|13/month  (;?36/y&ar).  Propane  is  piesently 
available  at  a cost  of  30.2'r’/gallon  ( 31  k'.Vh),  This  cost  will 
surely  rise  in  coming  years  causing  the  cost  of  the  back  up 
to  rise.  Therefore,  the  system  siiould  ba  designed  so  as  to 
iociuiro  as  little  back  up  as  possible  during  a "typical"  year. 

It  is  assumed  that  there  exists  an  exhausted  gas  well,  and  piping 
already,  thus  their  costs  are  not  considered, 
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Figure  4,1.2  - Graphical  Outpi 


than  the  allov/ed  amount  of  time  per  day,  the  back  up  (tractor 
enei'gy  in  this  case)  is  activated.  This  is  shown  between  the 
7th  and  8th  days  in  Figure  4.1.1.  This  use  of  tractor  energy 
costs  money,  which  is  reflected  in  the  increase  in  the  total 
capital  cost  of  tl.o  system. 

Figure  4.1.2  shows  the  aoouiaulation  of  "surplus"  energy, 
V^nenever  the  flywVieol  is  filled  to  capacity  with  energy  and  the 
wind  is  generating  additior.al  energy,  this  surplus  energy  is 
noted  by  TE3WI.  This  is  shown  between  liays  1 and  2 in  Figure 
4.1,2.  Whenever  the  wind  is  blowing  faster  than  the  rated  wind 
speed  of  the  system,  the  energy  available,  but  not  extracted  by 
the  system,  is  noted  by  TESW2  (see  Figuie  4.1.2). 

Figures  4.1.3  and  4.1.4  show  the  graphical  output  of  t\/o 
2 week  periods  in  Option  2,  Figure  4.1,3  shows  a period  from  the 
2r.d  quarter  of  the  year  (December,  January,  February).  During 
this  period,  no  itrigation  is  taking  place.  Therefore,  all  the 
energy  extracted  from  the  wind  is  going  Into  the  production  of 
hydrogen.  Whenever  the  wind  is  blowing  faster  than  the  rated 
wind  speed  of  the  system,  the  energy  available,  but  not  extracted 
ty  the  system,  is  noted  by  TESW2, 

Figure  4.1,4  shows  the  same  2 week  period  of  time  as  Figures 
4.1.1  and  4,1.2,  Irrigation  is  taking  place  during  this  period, 
therefore,  hydrogen  is  being  used  at  a steady  rate.  If  the 
hydrogen  storage  level  were  to  drop  below  the  low  energy  level, 
the  back  up  (propane)  would  he  activated.  Use  of  propane  causes 
an  increase  in  the  total  cost  of  the  system,  Si.nce  the  hydrogen 
storage  level  is  well  above  the  low  cnor.gy  level,  no  propane  is 
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U'jcd  during  the  w^ek  psriod  shown  in  Figure  4.1.,4. 

4.2  BASIC  EC0N0;-;iCS  OF  THE  SYSTEMS' 

The  final  total  capital  cost  at  the  end  of  a simulation  run 
contains  t'wo  costs,  the  cost  of  the  system  itself  and  the  cost  of 
operating  the  system  for  one  year.  The  cost  of  the  back  up 
systems  constitute  the  operating  cost  for  one  year,  while  the 
rest  of  the  total  capital  cost  is  the  cost  of  the  system  itself. 
Maintenance  costs  are  assumed  to  be  negligible  and  ai'e  not 
considered.  Taxes  are  also  not  considered  here.  These  cost  are 
suir.max'ized  in  Table  4.2,1. 

OPTION  1.  The  construction  of  the  system  modeled  by  Option  I 
requires  $145,000.  It  is  assumed  that  a farmer  ’,.'ishir.g  to  install 
Such  a systeiii  could  boi'row  the  necessary  capital  at  10>2  annual 
interest. 

Case  1.  In  case  1 the  lifetime  of  the  system  is  assumed  to  be 
10  years  with  zero  salvage  value  at  the  end  of  that  time.  It 
is  also  assumed  that  the  farmer  will  make  12  monthly  payments 
each  year  to  the  bank  against  his  loan  for  the  system  installation. 
The  amount  of  each  raonthly  payment  is  calculated  making  use  of 
the  formula  for  the  Capital  Recovery  Factor  (21)  given  in  KQN 
4.2.1. 

A = P EQN  4.2.1 

(1  + i)"  - 1 

where  A = end-of-period  cash  flo-ws  in  a scries  continuing  for 
a specified  number  of  periods; 

P = prcjjont  sum  of  money.  The  equivalent  worth  of  one 
or  more  cash  flo'ws  at  a relative  point  in  time 
called  the  presonti 
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i " r>o  (jffective  interoat  rate  per  period) 

N ==  -/..e  number  of  ccmpcunding  periods. 

V/ith  P ~ ':.i45,000 

i - 10,^  per  year/12  months  per  year, 

N = 20  years  * 12  months  per  year, 

then  A - $1 ,9l6/tnonth  = $22, 994/year 

The  total  cost  of  the  system  for  one  year,  then,  isj 

$ 22,994  - capital  recovery  cost 
+ $ 1.320  - operating  cost 

$ 24,314 


The  number  of  kilov;att-hours  of  energy  produced  during  this 
sane  j^erioJ  isi 

100  kV/  rating  * 90  days  ^ 24  hours/day  LOADPAC  = 131,^95  itA'h 
Dividing  the  total  cost  of  the  system  for  one  year  by  the  number 
of  kilov/att-hours  produced  in  one  year  yields* 

$24,314/131,695  l^Wh  = $0.lB/’ 

Casa  2.  In  case  2 the  lifetime  of  the  .3  is  assumed  to  bo 

20  years.  The  v/indmill  and  flywheel  system  should  actually 
have  a lifetime  of  more  than  this,  but  the  well  ai'^d  p\imp  may 
have  a lifetime  of  only  10  to  15  years. 

Using  EQN  4.2.1  yieldsi 

A = $1 ,399/:nonth  --  $l6,79l/yoar 
Therefore,  the  total  cost  of  the  system  for  one  year  isi 

$ 16,791  - capital  recovery  cost 
+ $ 1,320  - operating  cost 

$ 18,111  • • 
which  comes  to  a cost  of  $0.l4/kWh, 


OFI'ION  2.  The  construction  of  the  system  modeled  by  Option  2 
requires  $109,260.  Again,  10>^  annual  interest  is  assumed. 
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Case  1.  In  case  1 a 10  year  lifetime  v/ith  zero  salvoRC  value 
is  assumed,  Usitv'  EQN  yioldsj 

A - $1  ,444/inonth  - 4^1 7 1 327/year 
The  total  cost  of  the  oystiom  for  one  year  comos  toi 

$ 17,32?  - capital  recovery  cost 
f 26  - opera tir,'"  cost 

$ 17,353 

The  total  number  of  kilov.’at t--hours  produced  during  the  same 
period  of  time  iss 

60  k'lV  * 90  days  ’■  24  h.curs/clay  - 129,(^00  kV/h 
This  results  in  an  energy  coat  of  $0.1 3/k'.'/h. 

Case  2.  In  case  7 a 20  year  lifetime  v/ith  zero  salvage  value 

is  assumed.  Using  lil’J  4.2,1  yieldSi 

A - $1  , 034/r'.onth  - $12 ,658/year 

The  total  cost  of  the  system  for  one  year  oqualsi 

$ 12,656  - ojrital  recovery  cost 
t $ 26  - or:..ating  cost 

$ 12,684 

which  comes  to  an  energy  cost  of  $0.10/kiyh. 
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Cli.M'TKR  5 - Ut'SCL’oSlOr!  AND  COWMRNTS 
5.1  WiiY  Thc:3li  KOHELS? 

Only  two  sy3 tom  models  are  prosented  in  this  paper,  Tl’.ero 
are  a number  of  Gystom  models  t)iat  can  be  formulated  in  attempt- 
ing to  design  a system  that  would  iiarness  the  wind.  Option  1 
was  formulated  to  be  as  oinpl.e  a design  as  possible,  conver tiiig 
the  v;ind  power  to  mechanical  power  and  usirvi  this  mechanical 
power  directly  to  accomplish  tlie  task  of  pumping  irrigation 
wo  tor.  Starting  wimh  this  basic  system,  many  modifications 
are  possible.  The  flywheel  storage  could  be  eliminated  and 
replaced  wirh  a large  water  tank  or  reservoir,  'i’lie  system  could 
then  pump  water  year  I'ound  whenever  the  wind  was  blowing,  storin'? 
surplus  for  periods  of  low  wind.  Or  add  an  alternator  to  the 
system  and  you  can  drive  the  ii'rigation  pump  with  an  elo'itric 
motor.  During  the  non-summer  months  the  electricity  produced 
could  be  used  for  grain  dryin,5i  space  heating,  and  otlier  demands 
presented  on  the  farm.  Electric  utility  off  peak  power  could 
serve  as  hack  up.  Add  an  electrolysis  unit  to  the  system  and 
the  system  would  bo  essentially  that  presented  by  Option  2, 

Option  I was  chosen  as  a system  to  be  modeled  duo  to  its 
relative  simplicity.  The  balancing  of  the  rotor  size  and  fly- 
wheel (for  eiiergy  storage)  'size  deserves  a more  systeijatic 
api'roacn  than  is  presented  hero. 

Using  electricity  produced  by  the  wind  to  power  irrigation 
pumps  has  some  advantages.  However,  in  order  to  provide  the 


necessai'y  back  up,  the  irrigated  field  must  be  located  near  a 
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throe  plia.'jo  traroanit::: Ion  lino  to  invniniizo  tiio  coat  of  ili :: tribu tlon 
linen.  For  t!>e  naino  roanon  the  field  i;ina  t alno  be  located  close 
to  other  fo.rm  building,;:',  if  the  oleotricity  produced  durin,'  non- 
irrif.allon  periods  is  to  bo  uaed.  Tlic  pooaibility  of  uninj' 
exin  till"  power  company  networks  as  an  "infinite"  stora/;c  device 
dooa  exist.  This  would  probably  require  p/n’o rti.r.ent  intervention 
on  the  state  or  federal  level,  iiowever,  before  the  utility 
companies  and  other  parties  involved  could  r'acli  an  afjrecMont. 

Option  2 has  roveral  advantafjes  ever  the  otb.er  possibilities 
discussed.  The  r.iidwost,  where  it  is  proposed  to  install  such 
wind  powered  irrigation  systems,  is  blessed  'wit!'.  n ■ ny  natural 
pas  fields.  Natural  pas  is  presently  ured  to  power  ai'proxi:.ia  tely 
8'>,'  cf  the  irrigation  pui.ies  in  westorn  Kansas.  In  years  to  come, 
irore  and  mere  of  these  f;as  fields  '.vill  bo."in  to  riui  out  cf  natural 
pas.  There  are  already  numerous  depleted  wells  in  the  area. 
Fxistin"  pas  e.n.ines  and  pipelin'‘s  could  bo  used  In  Option  2. 
Depleted  wells  are  available  at  nominal  cost,  and  the  turbine 
could  bo  used  year  round.  By  expand  iiqj  the  system,  hydrojijcn  could 
be  produced  in  larco  cjiouqh  quantity  to  be  used  an  the  primary 
energy  source  on  the  farm.  The  idea  of  a liydropon  ecor.omy  is  not 

new.  (22,23). 

When  ex.iminir."  Option  2 it  should  be  pointed  out  tliat  this 
system  could  not  be  built  today  for  the  cost  cited  in  the  model. 
The  primary  obstacle  to  this  is  the  larpe  electrolysis  coll 
required.  The  specifications  used  In  this  paper  reflect  191^5 
projected  cells. 


5.2  WliY  TIIK'JE  NUMllKKS? 

It  has  already  been  pointed  out  that  the  selootion  of  the 
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i rotor  Dir.o  - i'^^fwhocl  r.izo  combiir'i lion  in  Option  1 licocivon  a 

; noro  nop'h i ■:  t icateii  approacli  tirxn  DYUAItO  can  tiivo,  Hov/ovor, 

• nevor.'l  di  'jioront  co.r.bini tLotvr.  woro  conaidorod.  None  of  those 

ooiahinatli  . ^ placoJ  “too  larfejo"  a (.■■oporrioncy  on  the  back  up 
system,  "Voo  lart;e‘'  is  fairly  arbitrary.  It  was  felt  that  the 

■ farxer  should  not  have  to  use  the  bad:  up  wore  than  once  or 

twice  a w('i;k.  From  Fi,3ui'e  5,2.1  it  can  be  soon  that  rotor  size 

< 

, of  1100  square  aetci's  witn  a flywi'.eel  enerpy  cap:icicy  of  1000 

ki Iowa tt-!'iCurs  is  about  in  th.e  middle  of  the  conibiraitions  that 

■ provide  a load  factor  of  approivima toly  6C;’>.  I t was,  tlieio fore , 

’ 2 

decided  to  use  a rotor  size  of  1100  m'  and  a flywheel  with  1000 

; kV/h  of  E torapio  capacity  In  tlie  simulation  of  Option  1. 

I 

; In  both  Option  1 and  Option  2,  the  cost  of  the  b.ick  up  is 

: actually  less  than  the  cost  of  operating  the  system  itself. 

This  enforced  the  desire  to  use  as  little  back  up  as  possible 
_ .1  'noneconomic*  dec  i si  on  must  bn  n id ' to  use  as  littli'  l-aok  up  as 

I possible  since  an  attempt  to  simply  minimize  tiio  cost  would  result  in 

! the  size  of  the  rotor  increases,  the  cost  of  tiie  system  increases 

linearly.  Therefore,  it  is  desired  to  have  t!io  rotor  as  small 
fis  possible  while  still  usin^  as  little  back  up  as  possible. 

Figure  $.2,J  shows  tliat  a rotor  size  of  700  square  meters 
. requires  very  little  propiano  bo  used  as  back  up  iii  Option  2, 

j Tliereforc , a rotor  size  of  700  square  meters  was  used  in  the 

J simulation  of  Option  2.  Figure  5»2.3  is  based  on  11)05  wind  data. 

Tins  data  was  deternuned  to  be  "typical’’  compared  to  averages 

( 

I over  25  years. 


5.3  WHY  A 3YSTKW  DYWAyiCS  AFFrOACH? 

System  liynamics  has  been  criticized  by  n-any  (2h)  on  the 
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f^round:!  tt’.at  Iho  aaLiuc.ptiotif:  used  an  tnany  nodola  have  not  been 
I'elovant  to  econonic  tl'.eory  and  that  the  abacnco  of  propoe  esti- 
mation techniques  result  in  considerable  abstraction  from  reality 
(25) • Du^,  as  has  been  stated  already,  the  models  presented 
here  bear  no  claim  of  being  "the"  correct  models.  They  merely 
represent  the  best  models  .-vailable  at  the  time  given  the  avail- 
able knowledge. 

It  has  been  said  that  iiacrensed  promotior;  ;idds  to  awareness 
of  a produce  (26).  Such  is  the  cane  with  systen  dynamics  and 
the  use  of  D’f?l\KO.  Originally  developed  as  a tool  for  use  in 
solving  industrial  management  prcbi.ems,  system  dynamics  have 
new  been  used  to  model  a wide  variety  of  topics  (4, 5,6, 7), 

The  author  feels,  tliat  by  using  system  dynamics  in  an 
area  tnat  m the  past  lias  been  ignorant  of  the  subject,  ho  can 
expose  its  existonce.  It  might  be  argued  that  DYIJAMO  docs  not 
present  the  b.st  approacii  to  attacking  the  problem  of  r.cdoling 
wind  powered  irrigation  systems.  3ut  it  is  an  approach  that  can 
bo  taken,  And  it  is  the  author's  hope  that  other  people  working 
in  these  areas  will  become  aware  of  the  potential  of  system 
dynamics  as  a method  for  problem  solving. 

5.4  A LCCK  TO  THE  FUTURE 

The  energy  costs  of  Option  I and  Option  2 compared  with  present 
energy  costs  are  shewn  in  Table  5*4. 1.  As  Ccin  be  scon,  even  the 
lowest  cost  associated  with  the  proposed  systems  is  2}  times  as 
high  as  the  higher  of  tie  two  present  energy  costs.  But,  Icok 
to  the  futxire. 

By  the  year  19()5  construction  of  Option  2 should  be  possible. 
Natural  gas  prices  will  have  surely  risen  considerable  by  then 
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Option  I 

1 Option  2 

Present 


Table  - Energy  (Worsts 

cost  (ip/\\‘lh) 


(10  year  loan)  O.lt) 

(20  year  loan)  0,14 

( 1 0 year  loan)  0,13 

(20  year  loan)  0,10 

electi'icity  0,04 

natural  gas  0,01 


due  to,  if  nothing  clno,  the  relative  ccarcLty  of  the  commodity. 
The  coat  of  fuel  to  utility  companiofj  will  likewise  have  risen 
rocultin,^  in  a hif^hor  cost  to  the  consumoi-  foe  electric  enorfiy. 

The  cost  of  construe tinjj  the  systems  modeled  will  similarly 
rise,  Eut  since  these  systems  use  the  free  wind  as  their  primary 
source  of  onor/jy,  the  rise  in  erei'cy  costs  for  Option  1 and  Option 
2 should  not  bo  as  great  as  the  rise  in  price  of  present  forms 
of  energy.  It  is  not  unlikely,  then,  that  v/ithin  the  next  10 
to  15  years  systems  such  as  those  px'uscnted  here  will  produce 
energy  at  a lower  cost  than  present  utilities.  Add  to  that  the 
fact  that  our  natural  resources  such  as  gas  arid  oil  will  not 
last  forever  , but  the  v/ind  should  continue  to  blow  always. 
Therefore,  it  would  not  be  unx'easonable  to  expect  these  types 
of  energy  producing  devices  to  become  wide  spread  in  the  future. 
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A IlHOKT  oumi-:  TO  DYNAMO 

1.  VAKtAUIK  NAMKS 

Varliblc  namco  ooruisl  of  fi'oiii  one  throvij'/n  rjcvoii  alphnbotic 
or  nunorlc  cbai'nc  tor'j , t);c  fLrr.t  of  which  trust  be  alph.ihij  tic , 
i'xainplo  I 

q 

INV 

LIST? 

7.,  KOUATION  TYTliS 

There  are  two  x-.iiii  variables  In  a dynaiiiio  model;  level 

variables  ar.i  rate  variables.  Aocorcii.r.^ly , tlie  oqvtations  Jofinittg 

these  variables  aro  called  level  and  rate  e.iuation3.  There  are 

other  .types  of  eiuationr. , too,  ' ach  as  auxilis.ry,  supplementary, 

and  ini  tial -value  of,u.ation3 . When  writin;'  a model,  the  type  of 

equation  should  be  spool  fled  in  the  first  column  of  Ihe  punched 

card.  Example  I 

card 

col 

I23'>i6?3. . . 

L LI  .K=I.l  .J+OTniNEL.JK-OTKL.JK) 

R RATT.KL=L:';V.K/nELT 

Other  equation  typos  are  spooified  as  followai 

A Auxiliary 

S Supplcmontary 

N Initial  Condition  or  Initial-Value 

Con-.-- tants . Constant.^  arc  di:sip;natcd  by  the  letter  C In 
the  first  column.  Several  constants  may  bo  specified  on  a card. 
Kxi\mpl  0 1 


/ 
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! card 

col 

12  3‘’»567«... 

i C DTA:^15 

I 

; C A--5.CAT-13.00.H1’1'=5.1251C-? 

I 

j P'^iiatlon  V.'r'l  t i rv^.  Tlic  ba::ic  rorir.at  of  lui  o>iuation  Isi 

qiiaiUity  name  I’xproraion 

I 

The  arithmetic  operations  used  to  from  at;  exprera lot;  aroi 

t 

+ Addition 
Snbtraotloi; 

1 * M;iltlpllcatlon 

j / Dlvlslot; 

! Kxpottentiatiot;  in  not  av'a  liable  an  at;  operation,  i>nt  power:;  and 

I 

roots  can  bo  calctilated  ur'ln.<;  the  f;inctioi;s  available  It;  DYNAMO. 

( 

Whet;  more  that;  ono  cper.;tloi;  appear:;  within  an  expression, 
multi  pi  Icatloi;  and  divi:;lot;  arc  computed  before  addltioi;  and 
i subtractioi;  with  the  o;'der  of  operatlot;  beitvt  from  left  to  right, 

’ Paren theses  may  be  u;;od  to  alter  ll;va?:o  normal  " p;'oocdei;ce.s''  ai;d 

j to  elimlt;atc  the  need  for  the  symbol  * for  re  pro  sent  1 1;'^  tho 

I 

, multlplicatlor;  of  two  iniantlties.  lixamploi 

A»-(l'-tC)'^(0+E)  is  equivalent  to  A'’(Htf.l)(D+E) 

Subscript:; . Level  variables  are  :;lways  s I tv’lo-subscr i pted , 
such  as I 

Ll.K  or  LEV3.J 

■ The  same  Is  Into  for  auxiliary  at;d  Gupplemei;tary  variables,  but 

their  slrtTle  subscript  Is  al;vay3  K, 

* 

Rato  variables  are  double-subscripted,  such  asi 
' RATT.JK  or  KLOWOUT.KL 


Table  A1..1  :5hould  bo  helpful  In  Idenll  fj  l tvi  the  correct 
subscripts  (11). 
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I Cat-')  !'’•  ' 'P/'.  Ivi’-KiLion  oi'  quanlity  Lypc  ia  to  ho  punched  j 

in  column  . .no  equation  chou'.d  start  in  column  ?,  lio  blanks 

' are  allcv/o>i  ..'•iiln  the  ciiuatiou  or  statom’nli  if  there  is  a«\y  i 

1 blank,  the  c-  nler  i(^nore:-:  wli.it  comes  afti'i*  tl.o  blank.  Oor.monts 

! may  then  bo  v.;itten  on  .statement  cards  provided  they  arc  ptoceeded 

by  at  least  one  blank,  should  the  m.iterial  not  fit  on  the  first 

1 I 

curd,  one  or  more  continuation  cards  nay  bo  used.  Tliesc  arc 

i 

I'unchod  with  an  X in  column  t,  and  the  mtorial  is  continued  starting'  j 

in  column  If  the  ri.'torial  is  continued  it  need  not  i~o  throujjh  j 

I 

colunn  72  of  the  ori«;inal  card,  hut  may  be  broken  after  a quantity 
name,  number,  or  arithcr.etic  op.’rator.  The  unused  portion  of  the 
, card  must  be  blank,  i.e.  no  comments  are  pernitted  within  the 

inferiration. 

Nuti^rj_eal  Cor ’^uta  tioiis . In  OY.'IAV.U  there  is  no  provision  for 
fixed-point  arithmetic.  All  c on vu tat  ions  aro  carried  out  in 
floatinii-point  form,  t’he  r.umorio'l  values,  however,  may  bo 
specified  as  fi  xo  l-po  int  v.aluos;  the  compiler  converts  them  into 

I 

' floati!'.;;-point  values.  Up  to  0 di';;its,  in  addition  to  sign  and 

; decimal  point,  may  be  used  to  express  numorioal  values.  For 

! 

• very  large  or  small  numbers  the  number  may  he  expros.'-cd  as  a 

number  multiplied  by  some  power  of  10  bj'  writing;  the  number 
followed  by  the  leitu'  H and  the  desired  power  of  ten.  h'xamploi 
■ 82  billion  can  be  written  as  82E+9  or  as  82K9  and 

1 ton- thousandth  can  bo  written  as 

t 

3.  DIHIMTION  CAKDS 

! nircolion  cards  arc  used  to  specify  the  leivtth  of  the  simulation 

i 

HIM,  time  Interval,  priiatiiv;  and  plottin,*;;  instructions,  and 
quantities  to  be  printed  or  plotted.  They  are  as  followsi 

i 

t 


115 


I.ionti  I'lcn  bion  Card.  This  is  t)ie  Tirst  cai'd  of  atiy  model, 
which  provides  a title  for  the  model.  It  begins  with  an  astorisk 
{*)  in  t)ie  first  column;  the  title  cannot  be  wore  than  40  char- 
acters in  Icagth  nor  can  ajiy  word  be  more  than  3 letters  long. 
Kxample i 

card 

col 

12345678. . . 

* WPI  - OPTION  1 

Run  Card.  Kach  run  and  rerun  is  assigned  a run  number,  by 
which  it  is  filed.  The  run  card  is  the  last  card  of  any  model, 
and  the  run  number  can  have  up  to  8 characters.  Example  1 

card 

col 

12345678. . . 

RUN  STD. 
or 

RUN  RS=-700 


SPEC  Card.  This  card  provides  values  for  the  following  four 
parameters  necessary  for  any  simulation; 

DT  the  interval  of  TIME  betv/een  TIME.J  and  TIiVPi.K 

LENGTH  the  value  of  TIKiE  when  the  run  is  to  be  terminated 

I’RTPKR  the  interval  of  TIKE  between  each  tabulation  of  the 
results 

PLTPER  the  interval  of  TIKE  between  each  plot  output  of 
the  results 

Those  four  parameters  may  be  defined  on  the  same  SPEC  card  or  as 
constants  on  separate  constant  cards.  Example; 


I 1*1111  ijiiTiiPi 
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card 

col 

12345678... 

3HiC  DT= . 1 25/LiiNC;TH=90/PR'"KiR-l/PLTfER= . 125 
or  as 

C DT^.125 

C IJiNGTH=90 

C PR'i’PKR=xl 

C PLTPKR=.t25 

PRiriT  Card.  This  card  specifies  what  values  are  to  be 

printed  and  in  v/hat  form.  In  any  table  of  printed  values,  TIME 

is  automatically  supplied  by  the  compiler.  Up  to  l4  quantities 

can  bo  printed  in  a table  output  in  addition  to  TIME.  A PRINT 

card  such  asi 

card 

col 

12345678. . . 

PR I NT  ABC , DEF , LEV , R ATT 
is  equivalent  to 

PRINT  ABC/DEP/Lli:V/RATT 
and  will  result  in  a printed  output  such  as« 

TIME  ABC  DEF  LEV  RATT 


But  a PRINT  statement  such  asj 

card 

col 

12345678... 

PRINT  ABC, DEF/LEV, RATT 

will  print  the  results  in  an  output  such  asi 


TIME  ABC  LEV 
DEF  RATT 


r 


1 

I 


11? 


In  printli'if',  DYNAMO  ^;ive3  only  up  to  5 3 v.‘;nif leant 
II'  values  are  very  small  or  very  bif^,  they  are  scaled  by  some 
power  of  10.  The  scalin';  factor  will  appear  in  the  title  of 
the  tabulated  output,  Tiio  scalin/^  factor  will  be  automatically 
provided  by  the  compiler,  ilov/ever,  the  user  may  specify  the 
sc'ilin;  factor  if  he  wishes.  For  example i 

card 

col 

12.-3h5678. . . 

PHIWT 

will  print  the  values  of  A after  dividin,;  them  by  10^\  and 
choosirvt  2 si  t.ii f icant  decimal  places  after  the  decimal  point. 

The  values  of  B are  divided  by  10^,  and  b-  decimal  places  are 
chosen  to  the  right  of  the  decimal  point. 

FhOT  Cards.  DYNAMO  has  an  automatic  plottir.g  feature  which 
enables  the  user  to  plot  up  to  10  quantities  on  a single  graph. 
The  scales  may  be  chosen  by  the  user  or  by  the  compiler.  All 
plots  are  verses  TIME.  Each  quantity  should  bo  assigned  a 
character  by  which  that  quantity  is  to  be  represented  on  the 
graph.  Examplej 

card 

col 

12345678. .. 

PLOT  LEV=L/RATT=$ 

will  plot  quantities  LEV  and  RATT,  represented  by  L and  $ 
respectively. 

Scales  may  be  specified  by  the  user,  or  else  they  will  be 
chosen  by  the  compiler.  It  is  possible  to  specify  only  the  upper 
or  lower  limit  of  the  scale,  leaving  it  up  to  the  compiler  to 
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specify  the  other  one.  The  unspecified  limit  is  given  as  as 

asterisk  (*).  Examplei 

card 

col 

12345678. . , 

PLOT  ABC=*/DEF=D( 0 , 200)/LEV=L( 0 ,* ) 
will  plot  variable  ABG  to  a scale  chosen  by  the  compiler  using 
the  * character,  DEF  to  a scale  from  0 to  200  using  the  D char- 
acter, and  LEV  to  a scale  with  a lower  limit  of  0 and  an  upper 
limit  to  be  chosen  by  the  DYNAMO  compiler  using  the  character  L. 

4.  COMMENT  CARDS 

Comments  may  be  made  on  NOTE  cards.  Example: 

card 

col 

12345678.. . 

NOTE  FWE  - FLYWHEEL  ENERGY  LEVEL  (KWH) 

Comments  may  also  be  added  to  a statement  card  following  at  least 

one  blank.  Example: 

card 

col 

12345678. .. 

N L1=0  INITIAL  CONDITIONS  FOR  LI 

5.  RERUNS 

A model  may  be  rerun  several  times,  for  different  parameter 
values.  Only  constants,  table  values,  and  direction  cards  may 
be  changed  in  a rerun.  Each  rerun  should  also  be  assigned  a 
run  number. 

6.  FUNCTIONS 

DYNAMO  has  a number  of  built-in  functions,  and  also  has 


1 
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provisions  to  accept  user-written  functions.  The  built-in 
functions  are  drsoribed  below. 

Common  functions.  Common  functions  defined  in  DYNAMO  arei 


EXP{A) 

A < 174 

Expontial 

L0GN( A) 

A > 0 

Natural  Logarithm 

SQRT(A) 

A > 0 

Square  Root 

sin(a; 

A < S23,000 

Sine  Function 

COS(A) 

A < 823,000 

Cosine  Function 

inc  and 

cosine  functions 

are  used  frequently  to  generate 

functions  of  TIME.  Raisin^t  numbers  to  a power  can  bo  computed 
as  I 

Y=A®  is  the  same  as  Y-^EXP(  i3’^L0GN(A) ) 

Random  Number  Generator.  DYNAMO  has  two  random  number 
generators j 

1)  NORMRN(MEAN,SDEV ) generates  random  numbers  normally 
distributed  with  mean  equal  to  MEAN,  and  standard  deviation  equal 
to  SDEV. 

2)  NOISE( ) generates  random  numbers  uniformly  distributed 
between  -0.5  and  +0.5.  Note  that  the  parentheses  are  necessary 
but  that  there  are  no  arguments. 

Third-Order  Delays.  DYNAMO  has  two  third-order  delay  func- 
tionsi 

1)  DELAY31 IN.DEL)  is  a material  delay  and 

2)  DLINF3( IN, DEL)  is  an  information  delay  wherei 

IN  = input  to  the  delay 
DEL  = magnitude  of  the  delay 


PULSE  Function 
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PULSE  ( , ER3'!' , I MTV  L) 

proti'ices  i jiubje  train  of  hej.j’ht  HCIiT,  with  width  of  DT.  The 
first  pulco  will,  appear  at  time  EKST,  and  thereafter  at  regular 
intervala  of  ler^^th  liITVL.  Neither  HGUT  nor  INTVL  need  to  constant. 

RAL P Function . 

RA.'-1P{SLP,STRT) 

is  equivalent  to  Figure  A2.1i 

RAMP=0  if  TIME  < STRT 

TIME 

RAt.P=  3LP»DT  if  TIME  > STRT 

STRT 

SAMPLE  Function. 

S AMPLE ( X , I NTV  L , ISAM ) 

aets  S.AMPI.E  equal  to  X at  rsaniple  tiir.es  separated  by  intervals  of 
length  INTVL,  and  holds  the  value  until  the  next  samplirig  time. 

ISAM  is  the  ir.itlal  value  of  SAMPLE. 

STEP  Function . 

STEP(iiGuT,STTM) 
is  equivalent  to  Figure  A2.2i 

STEP-0  if  TIME  < STTM 

STEP^^HGilT  if  TIME  > STTM 

Both  liPrliT  and  STTM  may  be  variables. 

?*AX  and  MIN  Functions . 

MAX(P,Q) 

setsi 


MAX=P 


if  p < q 
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i 

MIH(P,Q) 

so  1 

1 

M.IN=P 

if  P < Q 

1 

MIN^Q 

if  P ? J 

! • 

CLIP  Function. 

CL1P(P,Q,R,S) 

1 SO  ts  t 

\ 

\ 

CLIPr^P 

if  R < S 

/ 

1 

\ 

CLIP-^Q 

if  R > S 

/ 

f 

S'.v'ITCH  Function. 

■' 

* 

E’.v’lTClK  .P.Q.R) 

1 setsi 

1 

i 

SV/ITCH=P 

if  R = 0 

SWITCH=Q 

if  R / 0 

1 

1 

1 

TABLE  Function.  It 

may  be  desirable  to  express  the  values 

of  ono  variable  in  terms  of  the  values  of  anotlier  variable.  This 
was  tne  case  in  both  Option  1 and  Option  2 for  the  wind  speed 
with  respect  to  TIME.  The  TAEIJi  function  serves  this  purpose. 

The  form  of  a table  look-up  function  isi 

TABIfi  ( TNAME , X , XI NITL , XFI NAL , X I NCR  ) 

where  I 

TNAWE  - name  of  tlie  table 

X - independent  variable 
XINITL  - initial  value  of  rand  of  X 
XFINAL  - final  value  of  range  of  X 
XINCR  - increment  of  X 


J 


I 


Kxampi u I 

Siippo;  o Ih.il  Ihrt  follDwin"  v.iluoa  of  anxUl.-it'.y  varlablo  Z 
aro  ‘civon  with  r(';-.po<'t  to  iivlopeivlcMt  variable  Xi 


o'an  bo  exp re.areii  aai 


tZ-3456.'«. . . 

A z . K-  t'A  H 1 [ 'i’NA , x . k . i n , 6c , i o ) 

T T N ^ Tvs  i.  -■ ) X/'i  O/i-i  5/  6 0/6  5/1 2 0 

A’oto  litat  a V in  ijolijinji  1 in  used  to  denote  a 'I'AI'l.l'l  fiinotion. 
Mnenr  i n I or  po  la  t i on  ia  ur.ed  by  ilYNAMC  to  nj.T.putu  tlio  vaViioa 
that  are  not  opooificaliy  ,‘Tivon  in  the  table. 

In  TAHLl','  function;;,  X can  not  oxi'oeui  the  ispccified  ratv;e 
witaoiit  renorn  t i n,'.:  an  error.  'I'lie  TAblib  fuiictior.  ext.endn  tlio 
extreme  vahiea  of  I h(!  ciopondent  vtiri.",blo  If  the  independent 
varlalMo  oxoeods  the  specified  ran, 'to.  In  the  .abov'o  <?xainple, 

if  'l’Ai*iiI,  is  <ioed  in  place  of  'I’Aiil.i;,  will  have  a value  of  IZi) 
for  all  valnofi  of  X .^'.reator  (.hat;  6o  and  a value  of  55  I’or  all 
valuoo  of  X lena  that;  10, 

i'l'Li'i'TJJ l‘'tine  t i on.  Thir,  function  exponentially  r.mootlia  a 

(luaritity,  and  has  the  follow!  n,i;  forini 


IIMOU'I’IK  I N.UbL) 


where i 


IN  " input  to  I'o  iJii'.oo tin'll 
DPI,  umootliitv',  ooivitant  Of  liolay 
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This  Vs  only  at)  abbrnv  i a tud  fi;uulc  to  DYNAMO.  Poi*  inoi'tf 
Inrcfn-.atlon,  rofor  to  I'oforonco  (\l). 


